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I. LECTURE 1 
Introduction 

Let me start with a few words of explanation about the title of these lectures, "Quantum 
fields at finite temperature, from tera to nano Kelvin" . Tera is a Greek word which means 
10^^, while nano means 10"^ (this latter word is presumably very familiar to you because of 
the ^^nano technologies"). The Kelvin is the unit of temperature. In dealing with systems 
with so vastly different orders of magnitude it is useful to think in terms of the corresponding 
energy scales. Energies are conveniently measured in electron-volt. To convert Kelvin into 
electron-volt, recall that 300 K is about 1/40 eV or, if you wish, 1 eV is about 120000 
K. So 10^^ K translates approximately into 100 MeV. This is an energy scale typical of 
(high energy) nuclear processes. At the opposite end, the nano Kelvin corresponds to a 
subatomic energy scale. It is for instance realized in experiments with ultra cold atoms, 
where the phenomenon of Bose-Einstein condensation has been observed. The lectures will 
discuss theoretical techniques that are relevant to the calculation of the properties of matter 
in these two extreme energy regimes. The beauty of theoretical physics is that the same 
techniques are indeed capable to provide an adequate description of some important aspects 
of these vastly different systems. 

That quantum field theory appears as an essential tool in the description of hot and dense 
matter, composed of quarks and gluons, is a priori natural: the dynamics of quarks and 
gluons is governed by a quantum field theory, called Quantum Chromodynamics (QCD), of 
which I shall say more in today's lecture. What is perhaps more surprising is that field theory 
is also useful for understanding the behavior of cold atoms. Atoms are objects which we can 
study in isolation, and their dynamics obey non relativistic many-body quantum mechanics. 
However, collections of atoms can undergo collective, long wavelength oscillations. By long 
wavelength, I mean a wavelength much larger than the typical distance between the atoms. 
Such long wavelength oscillations necessarily involve collectively many atoms. And these 
collective excitations can be described by (classical) field theory. These long wavelength 
phenomena provide the connection between the two topics that I plan to discuss, hot and 
dense matter composed of quarks and gluons, and the Bose-Einstein condensation of weakly 
interacting cold atoms 
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A central theme of our discussion will be that, in both systems, the effect of the interaction 
can be large, although the strength of the interaction between the elementary constituents 
is small. At first sight, this looks like a paradox. But the clue to resolve this apparent puzzle 
has just been mentioned: collective, long wavelength phenomena involve many degrees of 
freedom, and the cooperation of these degrees of freedom compensates for the weakness of 
the coupling. Technically, this feature shows up in infrared divergences in the Feynman 
diagrams of perturbation theory. Because of these divergences, perturbation theory in fact 
breaks down, and other techniques have to be developed to perform calculations in the weak 
coupling regime. As we shall see, in both systems that we shall consider, a simple effective 
theory will allow us to overcome the difficulties met in perturbation theory (albeit only 
partially in the case of QCD). 

After this brief and general introduction, and the explanation of the title, let me say a 
few words about the plan of the lectures. There will be six lectures. The first four lectures 
will be mostly devoted to the physics of hot and dense matter, with in mind the quark gluon 
plasma. That will give me the opportunity to introduce techniques of quantum field theory 
at finite temperature which can be used in many other contexts. In these lectures, I shall be 
mainly concerned with the calculations of thermodynamical quantities (like the pressure). 
I shall use the scalar field as a prototype of a quantum field theory in order to illustrate 
the main difficulties that one encounters in perturbative calculations at finite temperature. 
Many of these difficulties are common to other field theories, in particular QCD. In the latter 
case further complications arise, that I shall briefly indicate. I shall also introduce a simple 
effective field theory that allows us to handle the infrared divergences of perturbation theory. 
This effective theory can be extended to the case of QCD, but this will not be covered in 
the lectures. We shall rather find a direct application of this effective theory in the study of 
Bose-Einstein condensation, to which the last two lectures will be devoted. Bose Einstein 
condensation is a phase transition which takes place in the ideal Bose gas, that is, without 
interaction. The question that I want to address is how the interactions between the atoms 
modify this remarkable phenomenon. More specifically, I shall be interested by the shift in 
the transition temperature caused by a very small repulsion between the atoms. Naively, 
since the interaction between the atoms can be chosen as small as desired, you may think 
of using perturbation theory in order to calculate the shift in T^. But we shall discover that 
perturbation theory is meaningless, it is plagued by infrared divergences. And we shall see 
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that the effective theory introduced in the first part of these lectures can be used to obtain 
an elegant solution to this problem. 

Finite temperature aspects of many-body physics or quantum field theory are presented 
in a number of textbooks, for instance [1-7]. Complements to the present lectures can 
be found in several lecture notes or review articles. Thus, the lectures 1-4 are based more 
specifically on [8, 9, 11-13]. Lectures 5-6 are based on [14]. These papers should be consulted 
for systematic references to the original literature. Further references will be given in the 
text about specific results that will be used or referred to. 

A. A brief introduction to QCD and its symmetries 

Let me now begin the discussion of hot and dense matter. At sufficiently high temper- 
ature and/or density, one expects nuclear matter -the matter that makes atomic nuclei- 
to turn into a plasma of quarks and gluons, whose interactions are described by Quantum 
Chromodynamics (QCD). I shall then remind you a few basic properties of QCD that are 
important to understand the bulk features of the quark-gluon plasma. 

1. Quantum Chromodynamics 

As 1 just said, QCD is the theory that governs the dynamics of quarks and gluons. Quarks 
are spin 1/2 fermions, that I shall represent by a field f{x). The index / refers to the so- 
called "flavor". There exists six such flavors, denoted u,d,s,c,t,b, and the corresponding 
quarks have different masses: the masses of u, d are very small, of order 2-4 MeV, that of 
the strange quark is in the hundred MeV range, that of the charm quark is of the order of a 
GeV, the same for the bottom, 4-5 GeV, while the top quark is much heavier, ~ 170 GeV. 
We shall be mostly concerned with matter that results from the "melting" of neutrons and 
protons, that is, matter made of up and down quarks. In addition to flavor, quarks carry 
another internal quantum number, color. There are A'c = 3 different colors for quarks. 

Gluons are somewhat similar to photons. They are the modes of a vector gauge field 
A^. As photons, they are massless bosons, with spin 1 (and two polarization states). In 
contrast to photons, which are electrically neutral, gluons carry color charge (gluons exist 
in N"^ — 1 = 8 colors), and interact directly among themselves. 
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The interaction between quarks and gluons is coded in the QCD lagrangian which takes 
the following form. 

^ = -^Fr^;, + ^fiijp - mf)mf, (1) 

where mj- is the mass of the quark with flavor f,p = '^^D^, with the covariant derivative 
= + ig^iJLi with the gauge field. This is a non-Abelian gauge field, i.e., a color 
matrix, = A^^t"", where t"" is a generator of the gauge group SU{3) (in the fundamental 
representation). The field strength tensor reads 

Fr = d>^A: - + gfatcA^.A'-,. (2) 

The first two terms are just those you would get in Quantum Electrodynamics (QED). In 
a non-Abelian gauge theory you have the additional term, quadratic term in A, which is 
responsible for the interactions among the gluons {fabc are the structure constants of the 
gauge group). 

Now the main thing that 1 want to do here is to remind you how the interactions among 
quarks and gluons, and among gluons themselves, can be red off the QCD lagrangian. Let's 
look first at the quarks. Their interactions with the gluon field is contained in the term 
^'-f^Afj^'^, and I shall represent this interaction by the first diagram in Fig. 1. The strength 
of the coupling is g. Let us now turn to the gluons. If one would ignore the last piece of 
the field strength tensor (2), i.e., set g = 0, then (-F)jv)^ would be quadratic in the gauge 
potential. And a lagrangian which is quadratic in the field describes only normal modes, or 
free particles. These modes are what we call the gluons. When g ^ 0, the gluons interact, 
and the interaction vertices can be obtained by analyzing the [Fj^^Y, which, aside from the 
part quadratic in the gauge potential, contains also the product of d^A"^ and gfabc^b^cy 
which generates a three gluon vertex, proportional to g and to the derivative of the field (so 
that the strength of the three-gluon interaction is proportional to the momentum of one of 
the gluons). We have also a four gluon vertex that is proportional to g^. These vertices are 
displayed in Fig. 1. 

2. Symmetries of QCD 

The symmetries of QCD play an important role in the analysis of the phase diagram of 
hot and dense matter. I shall briefly review here these symmetries. 
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FIG. 1: Feynman diagrams describing the elementary QCD interactions. 

There is of course the color symmetry, a local gauge symmetry. The QCD lagrangian 
is invariant under (local) color rotations of the quark field, accompanied by an appropriate 
transformation of the gauge potential: 

-> hA^h-^ - ^hd^h-^ ip ^ hip D^tp hD^^!, (3) 

with 

/,(x) = e*^''(^)*". (4) 

This gauge symmetry is actually the guiding principle behind the whole construction of the 
QCD lagrangian (1). 

Color symmetry is responsible for color confinement, the fact that color charges cannot be 
isolated. The quarks in Nature combine to form color singlet states, the hadrons. However, 
when matter is at extremely high temperature confinement "disappears", during the so 
called confinement / deconfinement transition. In these lectures, we shall be interested in the 
thermodynamics of this deconfined phase, called the quark-gluon plasma. 

When the mass of the quark is strictly zero, there is another symmetry, called chiral 
symmetry. Let me digress here on an issue that will appear at various occasions. The mass 
of the quark is something which is not determined by QCD, but by physics at a higher 
energy scale, by the electro-weak physics of the standard model. So the mass of the quark is 
something which is not "negotiable". However as a theorist, I can play games and consider 
a theory which is not quite QCD, but which differs from it only by the fact that = 0. If 
I allow myself to do that, I observe that the lagrangian possesses another symmetry, chiral 
symmetry, which corresponds to the independent rotation of the left and right quarks. It is 
a global symmetry. 
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What I did for the quark masses can be done for all the other parameters which are 
around. I have told you that = 3, and indeed in Nature there are three colors of quarks 
and eight colors of gluons. That is an experimental fact. But it is sometimes instructive 
to consider theories which look like QCD, but in which Nc can be varied. In particular 
you can sometimes obtain interesting insights by considering the limit where Nc is infinite. 
The theory corresponding to A'^c = oo has many features similar to those of the theory with 
Nc = 3. But letting Nc going to infinity allows you to do calculations which you cannot do 
if Nc is finite. We can also play with the number of flavors, etc. 

Returning to chiral symmetry, we note that it is spontaneously broken in the vacuum. 
This feature has many consequences, for low energy nuclear physics in particular: pions, 
the lightest hadrons, can be understood as (pesudo) Goldstone bosons. Chiral symmetry 
is also important in characterizing the bulk properties of hot and dense matter, as it is 
restored by thermal fluctuations at high temperature. Thus, the quark condensate (^^I/), 
an order parameter for chiral symmetry, is non vanishing in the vacuum, but vanishes at 
high temperature. 

One last symmetry which I want to mention, again without going into too much detail, is 
scale invariance. The QCD lagrangian (1) is left invariant in the rescaling of the coordinates 
X — Ax, where A is an arbitrary number, accompanied by a corresponding rescaling of the 
fields {A^ — 7- X^^A^). This symmetry (which holds for massless quarks) is easily verified by 
noticing that the gauge field has mass dimension 1, and that there are no dimensionful 
parameters in the lagrangian. In particular the coupling constant g is dimensionless. This is 
in fact a rather remarkable property because, as you know, the QCD lagrangian is supposed 
to allow you to calculate the mass of the proton for instance. So how are you going to get 
the mass of the proton in GeV, if there is nothing in the lagrangian that "knows" about 
GeV ? 

At this point, I need to remind you a few things about renormalization and the running 
coupling constant. When we go beyond the classical level, and calculate physical processes, 
we need to take into account the effects of the short wavelength fluctuations and that usually 
leads to divergent quantities. In order to control these divergences, you need to introduce 
some cutoff, and this is where an energy scale enters. In the case of QCD this energy scale 
is called Aqcd- It's value is of the order of 250 MeV. The key point now is that the original 
QCD coupling constant becomes a "running coupling constant", i.e., it depends on the scale 
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of the processes that one considers, or more precisely on the ratio of that scale to Kqcdi 
according to the (one-loop) formula 

This formula has a remarkable consequence, called asymptotic freedom: when /i is much 
bigger than Aqcd, g{lj) goes to zero. In the high temperature quark-gluon plasma, the 
typical energy scale is of the order of the temperature. Thus, when the temperature is large 
compare to Aqcd the coupling constant becomes small. This is essentially the argument 
that leads one to expect that the quark-gluon plasma is a weakly interacting system at high 
temperature. 

Now returning to the scale invariance, one notes that the the symmetry that exists at 
the level of the classical lagrangian is broken at the quantum level. This manifests itself 
in particular in a so-called "quantum anomaly", that can be measured by the trace of the 
energy momentum tensor. That trace, which should normally vanish for a system of massless 
particles (reflecting the scale invariance), is given (at one- loop) by 

K = ^MF^'F,,). (6) 

The function (]{g) is called the beta function. It describes the variation of the running 
coupling with the scale fi. It is given by 



At finite temperature, (and after subtraction of the vacuum contribution), one can write the 
trace of the energy momentum tensor as e — 3P, where e is the energy density and P the 
pressure . This quantity is non vanishing above the deconfinement transition, and (slowly) 
goes to zero with increasing temperature, in agreement with asymptotic freedom. 

I shall end here the discussion of the symmetries of QCD, and how these can help to 
characterize the bulk properties of dense matter. I shall return briefly to the phase diagram 
later today. At this point, the main message that I want to leave you with is that, because 
of asymptotic freedom, one expects matter at high temperature to be simple: a weakly 
interacting system of quarks and gluons. Since the interactions are weak it is natural to 
try and calculate their effects using perturbation theory. We shall learn in these lectures 
that the situation is in fact more subtle. But anyway, before calculating the effects of the 
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interactions, it is important to recall some well-known properties of the non interacting 
system. 

B. Thermodynamics of relativistic particles 

Let me then remind you about the thermodynamics of free relativistic particles. This 
will also offer us the opportunity of a short reminder of basics of statistical mechanics that 
will be useful later at various points in the lectures. 

1. Some reminders 

As you know, the statistical description of quantum systems involve the so-called density 
operator V, which, for systems in equilibrium at temperature T = is of the form 

where H is the hamiltonian, Q is the conserved charge which can be the electric charge, the 
baryon number, the strangeness, etc, and fi is the associated chemical potential. The fact 
that Q is a conserved charge means that H commutes with Q, [H, Q] = 0. It follows that 
the the eigenstates of the hamiltonian can be classified according to the eigenstates of Q. In 
other words one can write 

H = En Itpn) , Q \^Pn) = Qn \^n) • (9) 

We can also rewrite the density operator as 

V = J2\i^n)PnM, (10) 

n 

where p„ is the probability to find the system in the particular eigenstate {ipn) of the hamil- 
tonian. We have J2nPn = 1- 
The partition function 

Z = Tt e-^(^-'^Q) (11) 

is the central object of most calculations, since most thermodynamical functions can be 
obtained from Z. In particular, the thermodynamic potential reads 

n = -T\nZ. (12) 
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(I am using natural units where ks = 1. In other words, I am measuring the temperature 
in unit of energy or mass.) The thermodynamic potential is also 



where 



n = E-TS - nN, 



E={H) = ^Tr [He-f"^"-^^^) , 



(13) 



(14) 



and similarly is the expectation value of Q, N = (Q), which is calculated in the same 
way. S is the entropy. 

S = -ksTr VlnV = -kB^Pn'^npn- (15) 

n 

It is a positive quantity because p„ is a positive number smaller than 1. You can verify that 



E 



d\nZ 



dp 



N 



d\nZ 



9(/3/i) 



(16) 



The relation between the pressure and the thermodynamics potential is {V is the volume) 

(17) 



dV 



VL is the a function of temperature, chemical potential and volume. It is proportional to the 
volume 

fi(r,/i,l^) = ra;(/i,T), (18) 
from which one deduces that Vt = —PV. 



2. Free particles 

Now I would like to consider free particles. Let me write a typical free particle hamiltonian 
in second quantization 

p 

Here p represents the quantum numbers for one particle. It is typically the momentum, 
but I am not going to separate the momentum, the spin, the color, etc, and just use p as 
a generic symbol for the set of all the quantum numbers that are needed to characterize 
entirely the state of a single particle. I am assuming that the hamiltonian is diagonal in this 
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representation. And Ep is just the energy for a single particle in the state which is labeled 
by p. The operators and ap are creation and annihilation operator and these satisfy 
commutation or anticommutation relations, depending on whether the particles are bosons 
or fermions. That is, for bosons we have 



and for fermions 

Opc], + a^ttp = 1. (21) 
Now we want to calculate the partition function. 

Zo = Tr e-^-^° = Tr (^e~^^p"^"^''^) . (22) 

I have to tell you a bit more about what is this trace here. We deal with free particles. 
The partition function involves a sum over states that have arbitrary numbers of particles: 
the state with zero particle (the vacuum), the states with one particle, with two parti- 
cles, and so on. How do you characterize the state with n (identical) particles? You can 
look at the different states of one particle. Assume that the values of p are all discrete: 
Pq,Pi,P2, ■ ■ ■ iPki • • • is the list of all the possible states that a single particle can occupy. 
To characterize a state with a large number of identical particles, it is enough to tell what is 
the number of particles which occupy each single particle state p. The operator hp = a],ap 
is the operator which counts the number of particle in the state p. Now, I can easily rewrite 
Zq in terms of this operator. I get 

^0 = ( n ) = n ( I (23) 

{up} \ p / p y np J 

where J2{np} '^^^^ "configurations", that is all sets of possible numbers rip 

(eigenvalues of the operator hp). I have to distinguish two cases, the fermions and the 
bosons. For fermions, you cannot put more than one particle in a given state, so rip = 0, 1. 
The sum is then very easy: 

Zo = n (1 + e^^'-) (24) 
p 

For bosons, Up can be any integer, and we get 

= n (1 + + +■■■)= n 
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Quark-Gluon Plasma 




FIG. 2: The expected phase diagram of hot and dense hadronic matter in the plane {hb,T), where 
T is the temperature and hb the baryon chemical potential 

Knowing the partition function, you can calculate all thermodynamic quantities according 
to the formulae recalled above. In particular let me call fp = {"hp). This is usually referred 
to as the occupation number. In equilibrium, there is not an exact number of particles in 
each individual state, but each single particle state is occupied with some probability, and 
there are fluctuations. The occupation number can be obtained from the formula 

d\nZ 1 



fp 



(26) 



with — for bosons and + for fermions. 



C. The quark-hadron transition in the bag model. 

I now return, as promised, to the phase diagram of hot and dense matter, with a short 
digression on a simple model that mainly exploits the formulae that we have just recalled. 
Further details on this model may be found in [10]. 

The phase diagram of dense hadronic matter has the expected shape indicated in Fig. 2. 
There is a low density, low temperature region, corresponding to the world of ordinary 
hadrons, and a high density, high temperature region, where the dominant degrees of freedom 
are quarks and gluons. The precise determination of the transition line requires elaborate 
non perturbative techniques, such as those of lattice gauge theories. But one can get rough 
orders of magnitude for the transition temperature and density using a simple model dealing 
mostly with non-interacting particles. 

Let us first consider the transition in the case where /.i^ = 0. At low temperature this 
baryon free matter is composed of the lightest mesons, i.e. mostly the pions. At sufficiently 



13 



high temperature one should also take into account heavier mesons, but in the present 
discussion this is an inessential complication. We shall even make a further approximation 
by treating the pion as a massless particle. At very high temperature, we shall consider that 
hadronic matter is composed only of quarks and antiquarks (in equal numbers), and gluons, 
forming a quark-gluon plasma. In both the high temperature and the low temperature 
phases, interactions are neglected (except for the bag constant to be introduced below). 
The description of the transition will therefore be dominated by entropy considerations, i.e. 
by counting the degrees of freedom. 

The energy density e and the pressure P of a gas of massless pions are given by: 



where the factors 3 account for the 3 types of pions (tt^, tt^, 7r°). 

The energy density and pressure of the quark-gluon plasma are given by similar formulae: 

£ = 37 ■ — + 5, 
30 

P = 37 ■ — - P, (28) 
90 ^ ' 

where 37 = 2x8 + |x2x2x2x3is the effective number of degrees of freedom of gluons 

(8 colors, 2 spin states) and quarks (3 colors, 2 spins, 2 flavors, q and q). The quantity 

P, which is added to the energy density, and subtracted from the pressure, summarizes 

interaction effects which are responsible for a change in the vacuum structure between the 

low temperature and the high temperature phases. It was introduced first in the "bag 

model" of hadron structure as a restoring force needed to equilibrate the pressure generated 

by the kinetic energy of the quarks inside the bag. Roughly, the energy of the bag is 

P(P) = ^P^P + ^, (29) 

where C/R is the kinetic energy of massless quarks. Minimizing with respect to P, one finds 
that the energy at equilibrium is E (Rq) = 4PVo, where Vq = 47rPQ/3 is the equilibrium 
volume. For a proton with Pq ~ 1 GeV and Pq ~ 0.7 fm, one finds Po/Vq — 0.7 GeV/fm^, 
which corresponds to a "bag constant" P ^ 175 MeV/fm^, or B^/^ ^ 192 MeV. 

We can now compare the two phases as a function of the temperature. Fig. 3 shows how 
P varies as a function of T^. One sees that there exists a transition temperature 

T, = ( B'/' ^ 0.72 PV^ (30) 
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FIG. 3: The pressure of the massless pion gas compared to that of a quark-gluon plasma, showing 
the transition temperature T^. 

beyond which the quark-gluon plasma is thermodynamically favored (has largest pressure) 
compared to the pion gas. For 5^/^ ^ 200 MeV, ^ 150 MeV. 

The variation of the entropy density s = dP/dT as a function of the temperature is 
displayed in Fig. 4. Note that the bag constant B does not enter explicitly the expression of 
the entropy. However, B is involved in Fig. 4 indirectly, via the temperature Tc where the 
discontinuity As occurs. One verifies easily that the jump in entropy density As = Ae/Tc 
is directly proportional to the change in the number of active degrees of freedom when T 
crosses Tc. 

In order to extend these considerations to the case where /i^ 7^ 0, we note that the 
transition is taking place when the total pressure approximately vanishes, that is when the 
kinetic pressure of quarks and gluons approximately equilibrates the bag pressure. Taking 
as a criterion for the phase transition the condition P = 0, one replaces the value (30) for 
by the value {QO/SYtt^/^B^/^ ^ O.TOfi^/^, which is nearly identical to (30). We shall then 
assume that for any value of fiB and T, the phase transition occurs when P{fiB,T) = B, 
where B is the bag constant and P{fiBiT) is the kinetic pressure of quarks and gluons: 

nA.B,r) = 5I.^r + f(r^ + iJ), (31) 

The transition line is then given by P{fic,Tc) = B, and it has indeed the shape illustrated 
in Fig. 2. 
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FIG. 4: The entropy density. The jump As at the transition is proportional to the increase in the 
number of active degrees of freedom 

The model that we have just described reproduces some of the bulk features of the 
equation of state obtained through lattice gauge calculations. In particular, it exhibits the 
characteristic increase of the entropy density at the transition which corresponds to the 
emergence of a large number of new degrees of freedom associated with quarks and gluons. 
One should be cautious, however, and not attempt to draw too detailed conclusions about 
the nature of the phase transitions from such a simple model. In particular this model 
predicts (by construction!) a discontinuous transition; but this prediction should not be 
trusted [10]. 

D. Perturbative evaluation of the partition function 

I now begin the formal part of the lectures, and start introducing the techniques that will 
allow us to treat the effects of the interactions. 

The direct evaluation of the partition function of an interacting system is rarely possible 
as this requires a complete knowledge of the spectrum of the hamiltonian H . Various 
approximation schemes have therefore been devised to calculate Z. We briefly describe one 
of them, the perturbative expansion. We assume that we can split the hamiltonian into 
H = Hq + Hi with Hi <^ Hq, and that the spectrum of i^o is known Hq lip'^) = \ipn)- 
For convenience, we assume that the states {ip'^) are also eigenstates of the charge Q, which 
allows us to treat Q on the same footing as Hq. Thus, in the following, we shall assume that 
the term —fiQ is included in Hq. 
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1. The imaginary time formalism 

We define the following "evolution operator": 

U{t) = exp {-tH) = Uo{t)Uj{t), (32) 

where Uo{t) = exp{—HQT). The operator f/(r) is analogous to the familiar evolution op- 
erator of quantum mechanics, exp {—iHt). It differs from it solely by the replacement of 
the time t by — zr. Because of this analogy, we refer to r as the "imaginary time" (r is 
real!). It has no direct physical interpretation: its role is to properly keep track of the 
ordering of operators in the perturbative expansion (indeed in a "classical" approximation 
where the operators are allowed to commute the time dependence disappears). The operator 
Ui{t) = exp (riJo) 6xp {—tH) is called the interaction representation of U. We also define 
the interaction representation of the perturbation Hi. 

i7i(r) = e"^°i7ie-"^», (33) 

and similarly for other operators. It is easily verified that Ui{t) satisfies the following 
equation 

^ = -Hdr)U,(r), (34) 

with initial condition f//(0) = 1. By solving this equation one establishes the following 
important formula 

= e-/3Ho T exp ( - drHi (r) ] , (35) 



I. ^0 J 

where the symbol T implies an ordering of the operators on its right, from left to right in 
decreasing order of their time arguments. Using this formula one can rewrite Z in the form 

r./3 

(36) 



Z = Zo(Texp|-^ rfr/7i(r)|)c 



{O)o = Tr ( '-—O ] . (37) 



where, for any operator O, 

\ % 

The operator Hi is usually expressed in terms of creation and annihilation operators a and 
Then the calculation of Z reduces to that of the expectation values of time ordered 
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products of such operators. Whenever Hq is a quadratic function of these operators, Wick's 
theorem apphes, and all expectation values can be expressed in terms of single particle 
propagators. A diagrammatic expansion can be worked out following standard techniques. 
The partition function can be written a.s Z = Zq exp(rc) where Fc represents the sum of all 
connected diagrams. 

2. Free propagators 

The study of the free propagators will give us the opportunity to add a few remarks on 
the structure of perturbation theory at finite temperature. Let us consider a system with 
unperturbed hamiltonian: 

Ho = ^ eptt^ap, (38) 
p 

which commutes with the particle number operator Q = J2p CLpO-p- We define time dependent 
creation and annihilation operators in the interaction picture: 

at (r) = e^^^aje-"^" = e'-^a^ 

ap(r) = e^^^ttpe"^^" = e'^^^ap. (39) 

The last equalities follow simply from the equation of motion da^^\T)/dT = [Hq, a^'^^r)], 
and the commutation relations: 

[Ho,al] = epttl [Ho,ap] = -epap (40) 

which hold for bosons and fermions. The single particle propagator can then be obtained 
by a direct calculation: 

Gpin - T2) = (Tap(ri)at (r2))o 

= e-^-(^i-^^) [e{n - r2)(l ± Up) ± npd{T2 - n)] , (41) 

where: 

Up = {alap)o = ^^J^^ , (42) 

and the upper (lower) sign is for bosons (fermions). The fact that G is a function of ri — T2 
alone may be viewed as a consequence of the fact that Hq is independent of (imaginary) 
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time. Note that since < ri,r2 < /3, one has — /3 < ti — T2 < (3. One can verify on the 
expression (41) that, in this interval, Gp{T) is a periodic (boson) or antiperiodic (fermion) 
function of r: 

Gp(r - /3) = ±Gp(r) (0<r</3). (43) 

To show this, note the following useful relation: 

e^'p np = l±np. (44) 
Thanks to its periodicity (43), the propagator can be represented by a Fourier series 

^p(^) = ^E^"""^p(^^-)' (45) 
where the cj^'s are called the Matsubara frequencies: 

UJu = 2z/7r//3 bosons, , , 

' (46) 
OJ^ = (2z/ + l)7r//3 fermions. 

The inverse transform is given by 

1 

G{iuj,)= dTe^'^^G{T) = - — . (47) 

Using the property 

^(^) = ^E""""'^ -/3<r</3 (48) 

and Eqs. (45,47), it is easily seen that G{t) satisfies the differential equation 

{dr + Ho)G{T) = 5{T), (49) 

which may be also verified directly from Eq. (41). Alternatively, the single propagator 
at finite temperature may be obtained as the solution of this equation (49) with periodic 
(bosons) or antiperiodic (fermions) boundary conditions. 

Remark 1. The periodicity or antiperiodicity that we have uncovered on the explicit form 
of the unperturbed propagator is, in fact, a general property of the propagators of a many- 
body system in thermal equilibrium. It is a consequence of the commutation relations of 
the creation and annihilation operators and the cyclic invariance of the trace. I leave it to 
you as an exercise to establish this general property. 
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Remark 2. The statistical factor can be obtained from the relation = ±G(r = 0~ ,k). 
In the approximation in which the sum over Matsubara frequencies is limited to the term 
cOn = 0, one gets from (45): rifc ^ ^. Such an approximation corresponds to a "classical field" 
approximation valid when the occupation factors are large. This approximation, typically 
valid for long wavelength (small k), should not be confused with the classical limit reached 
when the thermal wavelength of the particles becomes small compared to their average 
separation distance. In this limit, the occupation of the single particle states becomes small, 
and the statistical factors can be approximated by their Boltzmann form: 
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II. LECTURE 2 

I briefly remind you where we stand. Last time I started to explain to you how to do 
calculations in field theory at finite temperature, and I told you that we were going to use 
two languages. One is the operator formalism that in many cases offers the most direct 
physical interpretation. This is a convenient formalism if you know the Hamiltonian. But 
(relativistic) quantum field theory is more often formulated in terms of a lagrangian rather 
than a hamiltonian, and the most appropriate formalism is then that of path integrals. I 
shall remind you how the calculations proceed in the two formalisms. 

Then, what 1 intend to do today is to apply these formalisms to the calculation of the 
thermodynamics of scalar fields. I would like to show you, on the simple example of the scalar 
field, how we go on calculating various Feynman diagrams, and alert you on the difficulties 
which emerge in such calculations. That will take us slowly towards the problems specific 
to QCD. 

A. Calculation of the partition function 

i. Operator formalism (reminder from last lecture) 

We established last time the basic formula for the calculation of the partition function, 
Z = TT{e~^"). We showed that it can be written as an expectation value of a time-ordered 
exponential: 

Z = Tre"'^-^ = Zo ^Texp |-^'^rfriJi(r)|^ , (51) 

where /3 = 1/T is the inverse temperature, and the expectation value is taken with the 
density operator associated with the free hamiltonian Hq. That is, for any operator O, we 
have 

(0)o = ^Tr[e-^^°0] . (52) 

Obviously, if you choose O to be the identity, this formula tells you that 1 = (1)q = 
■^tTe~^^° , which indicates that Zq = tre~^^" : Zq is the partition function corresponding to 
the hamiltonian Hq. I have written H = Hq + Hi, assuming that Hi ^ Hq, in some sense 
that I shall specify later on. I remind you that Hi[t) is the interaction representation of 



21 




FIG. 5: Examples of Feynman diagrams in (f)'^ theory. 

Hi, i.e., Hi{t) = e'^^"Hie~'^^° . The expression (51) suggests a way to calculate Z, as an 
expansion in powers of Hi. This is perturbation theory in the operator formalism. 

Now, what does it mean expanding in powers of Hi7 After expanding the exponential in 
Eq. (51), you have to take the expectation value. If Hq, as most often, is a quadratic form 
of creation and annihilation operators, for instance Hq = Xlp ^p^^pO-p^ then the calculation 
that you have to do can be expressed in terms of Feynman diagrams. The lowest order 
Feynman diagrams for the 0^ theory that we shall discuss soon are displayed in Fig. 5, with 
Hi entering the vertices. So the left diagram will be of first order in Hi, and the right one 
will be of second order in Hi. The lines which join two vertices at times ti and T2 are free 
propagators Gp{Ti — T2). Such a propagator, as we saw last time, has two pieces depending 
on whether Ti is bigger or smaller than T2: 

Giin - r^) = e-^-(^^-^^) [e{Ti - T2){1 ± n^) ± ^(rs - ri)np] , (53) 

where Up is the statistical factor (occupation number) for bosons or fermions, and the plus 
or minus sign corresponds to each of these two possible situations. Remember also that 
Gp(r) with r = Ti — r2, is a periodic function of r. Thus, it can be expanded as a Fourier 
series 

^p(^) = \ll e-'--^Gp{iuJn), Gp(za;„) = \ . (54) 
where Un = 2mT / (5 is a Matsubara frequency. 

2. Path integral formalism 

Now, I would like to move to the path integral formalism. The path integral formalism is 
based on the famous Feynman formula for probability amplitudes. Let me first discuss the 



22 



case of ordinary quantum mechanics. Consider a particle in one dimension, which is located 
at point Xi a the initial time t = 0. Its motion is described by a hamiltonian H, typically of 
the form H = /2m + V{x). At time t, the state of the particle is given by 

e-'^*|xi) . (55) 

The probability amplitude that at time t, the particle is located at point X2 is 

(xsle-^^^lxi) . (56) 

Now what Feynman told us is that this can be written as an integral over paths x{t) such 
that x(0) = Xi and x{t) = X2'- 

(X2|e-*^*|xi> = / V{x{t)) e'!oih-^-'-yi-))dt' ^ (57) 

Jx{ld)=xi 

where we recognize in the exponent the action integral. Note that the notation in Eq. (57) 
is symbolic: most of the paths involved in the sum are not smooth paths. But I'm assuming 
that most of you have some familiarity with this expression. 

What I shall do now is to use this expression (57) in order to calculate the partition 
function Z. I shall rely on the analogy that I already pointed out last time - the analogy 
between the evolution operator e~*^* and the density operator e~^^ of statistical mechanics. 
So I'm going to substitute it — )■ r, with r real. Then I shall calculate the matrix elements 
(x2|e~''^'^|xi) according to the rule we used in order to write the expression (57), and find a 
path integral from xi to X2, where now the path is a function of r. Watch out how things 
get modified in the action integral: it is replaced by r , idt is replaced by dr, and dt by idr, 
so that the kinetic energy acquires a minus sign. A minus sign exists also in front of the 
potential energy, and I pull this overall minus sign out of the integrand. Then I get 

(x2|e-^"|xi) = / V{x{t)) e-/o^(l--'+^(-))'^-' . (58) 

Jx{ld)=xi 

where now x = dx/dr. This is the formula that will be useful to calculate the partition 
function. Pay attention to the difference with Eq. (57): the factor i in front of the action in 
the real time path integral (57) has been replaced by a minus sign. And there is a plus sign 
in front of the potential instead of minus sign. 

Now I return to the partition function Z = Tre~~^^ , which can be calculated as 

Tre-^H = ! dx{x\e-f'"\x) . (59) 
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The matrix element is given by the path integral (58), with x{0) = x(/3) = x, and we shall 
write, symbolically. 



Tre 



x(0)=x(/3) 

where I have introduced the notation 5'^; for the "Euclidean action" 



Se 



(60) 



-mx^ + V{x) \ dr . 



(61) 



Eq. (60) is the path integral expression for the partition function. It involves a sum over 
paths x{t) that are periodic in imaginary time: x(0) = x{t). 

The extension of this formula to field theory is easy. The scalar field theory is specified 
by the lagrangian £ 



C = \{d,ct>f-\m'ct>'-V{ct>) 



(62) 



where I explicitly extracted the mass term from the potential V{(f)). I will very often consider 
a specific form for V^(0), the so-called "0^ field theory", where 

A 



V^(0) = ■ 



The notation (<9^0)^ is a shorthand for 

I now move to imaginary time, and change it into r. Then C is changed to 



C 



m'(t)' + y(0) 



(63) 



(64) 



(65) 



_2' ^' 2' ^' 2 
You see that the change of dt — )■ —idr produces the same minus sign here as in Eq. (60), 
and I can write the partition function of the scalar field as the following expression: 

Z= j V(P{t, x) e~^^ . (66) 

<l>{lS,x)=4>{0,x) 

Let me remind you that (p is playing here the role of the coordinate and is a function of time 
t and the three-vector x. I am summing over all field configurations which are periodic in 
the imaginary time r, i.e., 0(/3, x) = 0(0, x). The Euclidean action Se is an integral over r 
from to /3, 



dr d X 



1(9.0)2 + i(V0)2 + lmV + n0) 



(67) 
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This is the basic formula which we are going to use in specific calculations. 
Remarks 

In the operator formalism we deduced the periodicity of the Green's functions or the 
propagators from an explicit calculation in Fock space, calculating the time-dependence of 
the operators and observing the periodicity. 1 left it as an exercise to you to prove that this 
periodicity is in fact quite general and exists for the full Green's functions (i.e., not only 
for the Green's functions of the non interacting system). In the path integral approach, 
the periodicity emerges directly from the boundary condition in Eq. (66): The trace in the 
partition function involves summing over field configurations that are periodic in imaginary 
time. 

There are further remarks which I want to make at this stage. 

• The first remark concerns the Euclidean metric, and why the action is called Euclidean. 
By Euclidean metric I mean that all the gradient terms come with positive signs. 
This makes the action Se a positive definite quantity, provided of course that the 
potential V{(j)) is well-behaved. For example, if it is A0^, we require A > 0. Then, one 
can interpret this exponential e"*^^ as a Boltzmann weight. That is to say, one can 
interpret e~'^^ as a probability distribution, which allows us, in particular, to calculate 
the path integral using Monte-Carlo techniques: One discretizes the field </>, puts it 
on the four-dimensional lattice and select the paths using "importance sampling", i.e., 
with a probability distribution proportional to e~'^^. You would not be able to do that 
in real time, because then, you would have to add wildly oscillating factors e*'^, which 
no computer knows how to do. 

• The second remark concerns the periodicity condition in the expression (66). Let's 
assume for a minute that we forget about it. What do we get? We may rewrite the 
Euclidean action as 



where V represents the gradient in four dimensions. This may be interpreted as the 
energy of a classical field configuration in four dimension, and the whole path integral 
can be viewed as the partition function for a classical field theory in four dimension 
(a sum over all classical field configuration weighted by the factor e""^^). So, if one 




(68) 



25 



ignores the temporal periodicity, one is left with a four dimensional statistical field 
theory. 

Q: In this case, isn't there a temperature factor in front of the action? 

A: Indeed, strictly speaking, the Boltzmann factor is the exponential of 1/T times the 
energy of the field configuration (here Se in Eq. (68)). However, this factor 1/T 
is just a multiplicative normalization. This is not the main point that I want to 
emphasize here, which is that finite temperature effects in a quantum field theory 
could be viewed as finite size effects in a problem of classical statistical mechanics 
in one dimension more. In fact, by abandoning the periodicity condition and 
letting the integration over imaginary time extend to infinity, one is looking at 
the zero temperature limit, that is, one is doing quantum field theory for the 
vacuum. 

• The third remark concerns again the integration over the imaginary time in Se- I 
have argued that, if I let /3 — )■ oo, I can treat the time r as an ordinary coordinate, 
and I end up with a four- dimensional classical field theory. Now, I want to argue the 
other way. Let's consider the limit where /3 — ?■ or T — t- oo, that is, the limit of very 
large temperature, large compared to all typical energy scales in the problem. For 
instance, when you have a mass m, this limit applies when T ^ m. If /3 is very small, 
and unless extremely singular field configurations play a role, I can ignore the time 
dependence of the field (p. Then I can rewrite the Euclidean action as 



Because the fields are considered to be independent of r, the r-integration can be 
done trivially and leads to the factor /3 in front. Now we get again a classical field 
theory, this time in three dimensions. The reason why physics is becoming classical 
here is because I am ignoring the time dependence. The role of the imaginary-time 
dependence is clear in the operator formalism. The time dependence comes because 
the operators do not commute. You can look at the equation of motion 




(69) 



dO 
dr 



[O.H]. 



(70) 
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If the hamiltonian commutes with the operator O, then O is time independent. That 
is, the imaginary-time dependence of an operator is related to the non-commutation 
of the operator in question with the hamiltonian. Thus, in the formula (51), the time- 
ordered exponential is there to keep track of the time-ordering of quantum mechanical 
operators. Now, I would like to remind you of one thing. Recall the formula for the 
propagator that I wrote at the beginning of the lecture 

Gl{r) = e--(— ) [e{T){l ± n^) ± e{-T)np] , (71) 

where t = ti — T2- Clearly, in the limit where r — )■ 0~: 

Gl{T ^ 0-) = ±n^ . (72) 

Let me focus on bosons for which 



1 

n. 



up 



p/T - 1 



e' 



(73) 



When ep/T is a small number, I can expand the exponential to obtain 



~ - . (74) 



ep 



This is the occupation factor that you get in classical field theory. 
Now remember also that G'p(r) can be expanded in a Fourier series 

The approximate statistical factor in Eq. (74) is obtained by ignoring all Matsubara 
frequencies except the one with n = 0. 

These are features that we are going to meet repeatedly. 

Q: You are assuming h = 1 here, but if you put ^'s back in this formula one of ^'s appears 
in front of 5*^, and I guess h is something like temperature... 

A: I am going to leave as an exercise to you to put back the h everywhere and understand 
how the classical aspects enter as /i — t- 0. But the classical physics which we are 
discussing is actually subtle. One sometimes joke, that there are "two" Ws: 

there is the h of Mr. Pauli and that of Mr. Heisenberg. These are of course the same 
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H, that of Mr. Planck (and we set it equal to one !). What I want to emphasize here 
is that there are two physics issues. We are going to deal with systems of particles 
which are intrinsically quantum, like the black body radiation. There the h involved 
is that of Pauli, that of the "Pauli principle" . And we are going to deal also with long 
wavelength excitations. These behave like classical field oscillations whose wavelength 
is large compared with the typical distance between the particles. That allows us to use 
the small gradient expansion and that expansion is controlled by the h of Heisenberg 
(that of the "uncertainty principle"). These features are often mixed in a subtle way. 
You can have long wavelength, classical-like, oscillations of a gas of particles which are 
intrinsically quantum. So, taking the limit — )■ is tricky here. 



B. Thermodynamics of the scalar field 

With all this preparation, we can now begin the discussion of the thermodynamics of the 
scalar field. I shall do that using perturbation theory, and shall use a mixed formalism, just 
to train you to go from one language to the other. Let me view the scalar field in terms of 
the hamiltonian to start with. This is 

H = ld'x Qvr^ + 1(V0)^ + + ^0^) = H, + H, , (76) 

where 7i{x) is the canonical momentum conjugate to 0, and ifi ~ |y J 0^. When I said 
earlier that Hi is small compared to i^o; it means in the present context that the coupling 
constant A is a small number. And I want to calculate the thermodynamic potential as a 
power series in A, 

n = -^\ogZ = no + \Qi + x^n2 + --- ■ (77) 

The reason why 1 write the hamiltonian first is that I don't want to spend time on 
Qq, which is associated with the free hamiltonian Hq. 1 can write Hq in terms of the 
normal modes of the field, that collection of harmonic oscillators with frequencies 

Uk = \/k^ + m? (corresponding to the wave numbers k): 

Ho = ^Uk (alak + ^ ) • (78) 
fc ^ ^ 
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FIG. 6: First order correction to the thermodynamic potential. 
As you probably all know, from this expression one can easily calculate Qq as 

k ^ k 

that is, the sum over all the modes of the corresponding oscillator thermodynamic potentials. 

1. Short wavelength modes 

This expression (79) reveals a problem which we are going to face repeatedly: ultraviolet 
divergences. I want to spend a few minutes on this. 

The sum over k for the zero point energy is a shorthand for 

As it stands, this integral has no meaning, because Vfc^ + ~ |fc| for very large k and if I 
put a cutoff A at the upper end of this integral the result will grow like A^. So it is an infinite 
number when A — )■ oo. On the other hand, it does not depend on the temperature T, and can 
be interpreted as a correction to the vacuum energy. Accordingly, we are going to redefine 
the zero of the energy and simply subtract it. This is a "poor man's renormalization" and 
we will do more sophisticated things soon. What I am doing here is simply redefining the 
zero of the energy in such a way that the contribution of the zero point energies (the first 
term in Eq. (79)) drops out. The second term in Eq. (79) is finite because when -> oo the 
factor e~^'^'' goes to zero. 

Now I want to move on and proceed with the calculation of the first-order correction, Qi. 
It is given by the diagram in Fig. 6. Because this is the first order of perturbation theory, 
this is simply A/4! times the expectation value of 0^ calculated in the non-interacting ground 
state. Since the non-interacting energy is quadratic in the fields, the average is a Gaussian 
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integral, and therefore (0'*)o = 3 (0^)q (0^)o- Now, (0^)o is nothing but the propagator 
evaluated at r = and x = 0. Thus 



4! 



A 



>o = 7;3(0^>o = 9[^(- = O'^ = O)]'- 



A 



4! 



The boson propagator is, in a mixed representation 

1 



(81) 



AJr) 



2uji 



[{l+rik] 



] ■ 



(82) 



You can also verify that this is a periodic function of r and it has the Fourier decomposition 
like the propagators that 1 was discussing earlier. The Matsubara representation is 

1 



(83) 



The propagators I used before were those for the creation and annihilation operators, and 
they had the singularity at r = because it matters whether you order on the left of a 
or a on the left of a''". For the propagator of the scalar field, there is no singularity at r = 0. 
In particular for r = and x = 0, 

d^k 1 + 2nk 



A(r = 0,x = 0) 



(27r)3 2iUk 



(84) 



This is an important formula that we are going to use again and again. This is the expression 
of the fluctuation of the field, (0^)q in the absence of interactions. 

The calculation can be completed now. The expression (81) is equal to 



A 



1 + 2nfc 
2wfc 



^2uJ ^(?a;J ^^\^2uJ i^Uk 



(85) 



where I rewrite the integrals as sums over k. Actually there are various problems in this 
expression. 

Look at the first term: 

2 / „ \ 2 

(86) 



^ 2uJk 



1 



(27r)3 Jk^ + 



This integral is divergent. But since it doesn't depend on temperature, we can drop it, as 
we did earlier: it represents an infinite correction, of order 0(A), to the zero of the energy. 
The second term is finite, because of the presence of the statistical factor. 
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FIG. 7: First order correction to the mass. 



The troublesome term is the third one. This is divergent and temperature dependent. 
We know from general principles that this should not happen. What I want to show you is 
that it will disappear. To that aim, it is useful to calculate the leading order correction to 
the mass. This is given by the simple diagram in Fig. 7. Let me call that S. It is a number 
which doesn't depend on the momentum of the external lines and is equal to 



The propagator G^, including the mass correction, is now given by the Dyson equation 



where is given in Eq. (82). 

Now, we observe that, in S, there is a contribution which is independent of temperature, 
so which will exist in the limit of vanishing temperature. This contribution is divergent, but 
I cannot handle this divergence the way I did for the energy - just subtract it - because the 
correction to the mass propagates in all the lines of the Feynman diagrams. So I need to do 
a more elaborate treatment. I introduce a term, called a "counter-term". 



and add this to the hamiltonian density. Then I require that this correction doesn't change 
the mass at zero temperature, that is. 




(87) 




(88) 




(89) 




(90) 



This determines 6m'^: 




(91) 
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FIG. 8: Three-loop contribution to the thermodynamic potential. 



Of course, the term added to H contributes to the thermodynamic potential. It will 
in particular generate a contribution to Qi. Since 6m'^ is already of order A, in order to 
obtain this contribution, one just needs to calculate the expectation value of 0^ in the 
non-interacting system, which we know already (see Eq. (84)). Thus: 



2 ' ' V / 2 \ 2 ^ 2a;fe / \ ^ 2ujk 

\ k / \ k 

Now if you compare what we have here with Eq. (85), you have A/4 x (divergent sum) x 
(finite sum) there, and now A/4 with minus sign x (divergent sum) x (finite sum) here. 
So they cancel, as anticipated. This simple example illustrates a general result: there 
cannot be quantities which are (ultraviolet) divergent and which depend on temperature. 
Such contributions are eliminated by a proper treatment of the sub divergences of Feynman 
diagrams, of which I have given you a very elementary example. This property can be 
verified to all orders. The final result for Qi reads 

fii = -^T\ (93) 
1152 ^ ^ 

I would like to show you now the results of the three loop calculation, which involves in 
particular the diagram displayed in Fig. 8, in order to illustrate other important features of 
field theoretical calculations at finite temperature. When we do such higher loop calcula- 
tions, we need more sophisticated techniques for the renormalization procedure than what 
I've done so far. In particular in gauge theory it is essential, and in the scalar case it is 
convenient, to use dimensional regularization. The outcome is that, in higher orders, quan- 
tities like the coupling constant, or the mass, start to be dependent on the "renormalization 
scale" . This dependence is what I would like to discuss now. 
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FIG. 9: One-loop correction to the coupling constant A. 



Let me just quote the result for m = (fi = —PV): 



^ vr^r^ I 1 1 A 



9 10 8 167r2 8 



31og J^ + ^ + C 
^ 47iT 15 




(94) 



You could have guessed (from last lecture) the first term 7r^T^/90, which is just the pressure 
of non-interacting massless modes. C is a number which has an explicit expression in terms 
of dilogarithms and the Riemann zeta function, but this expression doesn't matter here. In 
this expression, fi is the renormalization scale which can be chosen at will; it can be 1 MeV, 
can be 1 GeV, can be whatever you want, a priori. 

Now, the pressure is a physical quantity. It cannot depend on what you choose for /i, and 
therefore dP/dfj, has to be zero: 

<-) 

How is that possible since there is an explicit fi dependence which just comes from the 
calculation of the diagrams? The point is that A depends also on /x: 

/.|-A(/.) = /3(A) = 3^ + 0{X') . (96) 

(Note that the /3-function is here positive, in contrast to QCD.) This dependence on n 
appears essentially because A receives contribution from diagrams, like that in Fig. 9, that 
are logarithmically divergent. The factor 3 comes from the three independent channels. 

Now let us calculate dP/dfi. This contains two contributions, the /i dependence implicit 
in A, and the explicit one. We have 

dP n^T^ f 1 1 dX 3 / A \'l 

Here, in line with the weak-coupling expansion, I am ignoring the /i-dependence of the A^ 
term, as this is of order A'^. By using Eq. (96), one verifies that the term within the braces is 
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zero, In other words, if you do a calculation at the order of three loops, which is a calculation 
at the order of A^, then the pressure P is independent of the renormalization scale /i at this 
order. That can also be verified in higher order calculation. It is a general result. 

But there is more to be said. You see that /i enters in a logarithm. This is also a fairly 
generic situation when doing high order calculations. Now, in the logarithm, /i is divided 
by some scale, and the typical scale that appears is that in Eq. (94), namely /i appears 
generically as yu/(27rT). The logarithm can be large if /i is very big compared to 27rT, or if 
/i is very small compared to 2'kT. You don't want that because you are doing an expansion 
in powers of A, and in order for the successive terms to be as small as possible, you want 
these logs to be as small as possible. This is why it is natural to choose ii of the order of 
2tiT. 

This is an important observation, which leads one to expect the thermodynamics of QCD 
at high temperature to be very close to that of a free gas of quarks and gluons. The reason 
is that, in QCD, there is a minus sign in the /3-function, so that the coupling decreases as ix 
increases. Since the optimized /i is related to the temperature as we have just indicated, at 
very large temperature /i is large and the coupling is small. 

This is the summary of the arguments which deal with short wavelength fluctuations - 
the modes with very large momenta. There are ultraviolet divergences that we can control 
by the standard process of renormalization. Once this is properly done, we get finite results 
at any finite temperature. I have shown to you one particular example where indeed terms, 
that may appear in intermediate stages of a calculation, which are divergent and depend 
on the temperature, do cancel. And finally I've shown that in higher order, you have to 
worry about running coupling constant effects. The running of the coupling is in particular 
essential to guarantee that the physical observables are independent of the renormalization 
scale /i, at the order at which we calculate. I have also indicated that a natural scale may 
be chosen, /i ~ 27rT, in order to optimize the apparent convergence of perturbation theory. 

2. Long wavelength modes 

I would like to go now into another regime, that of the long wavelength modes. I'm going 
to address another technical issue which is related to infrared divergences. This is actually 
where the major difficulties of perturbation theory at finite temperature lie, and this will 
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FIG. 10: A second order correction to the mass. 



also occupy us next time. 

To introduce the subject, let me return to my calculation of the correction to the mass, 
Eq. (87). I am going to drop systematically the vacuum contribution. I shall work in the 
limit m — where I can do explicitly the relevant integral: 



Note that this integral would be quadratically divergent in the absence of the statistical 
factor. (This factor provides a cut-off at the scale /c ~ T, and dimensional analysis shows 
that the integral is proportional to T^. ) Thus it is dominated by "hard" modes, with 
k = 0{T). The modes which contribute to the integral are those plasma particles which 
have wavelength of the order 1/T. This mass correction is actually important, and I'm going 
to call it the Debye mass m^. 



The inverse oirriD plays the same role as the Debye screening length in an ordinary plasma. 
This is an example of what we will meet later and what is called a hard thermal loop. This 
terminology has its origin precisely in the fact that the integral is dominated by hard modes. 
This correction is also sometimes called "thermal mass" . 

Let me now imagine calculating the second-order correction to the mass, given by the 
diagram in Fig. 10. What do I get? There is a minus sign because there are two vertices. 
Then I have a sum over the Matsubara frequency, and an integral over the momentum. 
The insertion is m|,. We get (I am assuming again that the particle is massless at zero 



2 J (27r)3^ ~2 J (27r)3e^-/r_ 1 A? ~ 2 12 ' 



(98) 
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FIG. 11: Mickey Mouse diagram for the mass correction. 





FIG. 12: Ring diagrams for mass correction (left) and thermodynamic potential (right). 



temperature): 
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2' (27r)3 (a;2 + P)2 " 

Now look at what happens in the particular term with n = in this sum, 
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(101) 



(27r)3 (P)2 • 

That is awful. Because k'^dk/k^ = dk/k'^, this is infrared divergent as J dk/k'^. We have a 
problem. 

And this problem is going to be worse and worse as we try and calculate more diagrams. 
Let's indeed calculate the "Mickey Mouse" diagram in Fig. 11. This is proportional to the 
integral 

d^k (mlY r dk 



(102) 



(27r)3 (a;2 + p)3 J (^2)2 • 

That's indeed worse! The more mass insertions you add, the more serious the problem is. 

But now something happens. You see that what we are doing here is something a little 
bit stupid, once you understand the physics. I will have the opportunity to come back to 
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this later. What we have recognized is that, in the thermal bath, the particle acquires a 
mass even if they are massless to start with. This is a non-perturbative effect, although the 
calculation of the mass itself can be done within perturbation theory. But the fact that the 
particle acquires a mass is something that you should keep in mind when you do a higher 
order calculation. When we do this calculation order-by-order, we extract the mass and we 
treat it as a correction. We should not do that. If the particle has a mass, let's take this 
properly into account! Let me do that by calculating the sum of an infinite number of such 
insertions in Fig. 12, which 1 will call ring diagram. You can write the expression for that 
because that is the same as Fig. 10 but with a slight modification: 

where the first propagator, 1 / (u;^ + fc^+m|,), is that of a massive particle with mass m^- You 
can verify that if you expand this propagator in powers of m|,, then you generate back all the 
diagrams that we have considered before and that are infrared divergent. The second term in 
Eq. (103) subtracts the contribution of the diagram in Fig. 7, which is calculated differently 
(there is no infrared divergence in this hard thermal loop contribution). The remaining 
integrand is then dominated by soft momenta. Now, for the particular contribution with 
n = 0, 1 get 

You see that the mass in the denominator provides an infrared cut-off that allows you to 
calculate explicitly the integral. Simple dimensional analysis reveals that it is proportional 
to m/) ~ -\/Ar. An explicit calculation, using the explicit expression for m^) given above, 
yields 

ATV A ^ 



Likewise, you can calculate the ring contribution, Qfmg, in Fig- 12 to the thermodynamic 
potential and you get 

I rj.\^ f d^k ( \ ml ] ml \ 



I rr. f d^k . 



1 + ^ 
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This is an important result. It tells you that if you do an infinite resummation, you get rid 
of the infrared divergences. Infrared divergences appear then an artifact of the fact that we 
did not do the calculation in a proper way, i.e., we attempted an expansion that is not valid. 
If we include, correctly, the thermal mass in the propagator, we get a finite result. 

The second lesson is that you see something very funny appearing. I told you before that 
we assume a power series in A: Q = Qq + XQi + X'^Q2 + ' ■ ■ • And this is indeed what you expect 
if you calculate Q by just calculating Feynman diagrams one after the other. But what you 
get here is a non integer power, X^^'^l An unexpected term is sneaking in the series, which is 
A^/^fi3/2. That means that the weak coupling expansion of the thermodynamic potential is 
not what you naively expected. You get a contribution which is not analytic in the coupling 
constant. And there is no way to get such a contribution by summing a finite number of 
terms. 

There is a lot of material in this lesson. We will come back on several points next time. 
I will then try to explain what is going on in very simple terms. 

Q: In this result for the thermal mass, A should be positive, shouldn't it? 

A: The A should be positive in any case, because the potential should be bounded from 
below for the system to be stable. 

Q: Is the instability related to the appearance of the A^^^ term? 

A: No, but indeed, when A is negative, you have to pay attention to the meaning of A^/^. 
But this is not really related to this instability issue that I am referring to. If the 
system is truly unstable, worse things happen. 

Q: If you have an attractive interaction then the mass term will be eaten up by the inter- 
action, which can give a vanishing mass. I guess that is what he is asking. 

A: The potential of the 0^ theory is like this (Fig. 13, right). It could also be like this 
(left of Fig. 13), corresponding to < 0. The effect of thermal fluctuations would 
be typically to restore the symmetry, i.e., to transform the potential from its initial 
shape (at T = 0) on the right hand side of Fig. 13 to the convex shape of the l.h.s of 
Fig. 13. That is, thermal fluctuations contribute positively to the mass squared, and 
if < to start with (at T = 0), the thermal fluctuations will eventually turn it 
positive at sufficiently high temperature. 
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FIG. 13: The potential of the cp^ theory wihtout {m? > 0) and with {m? < 0) spontaneous 
symmetry breaking. 

Q: Is mjj proportional to A? 

A: Yes, m'jj is proportional to A (at least its leading contribution), and it is positive if A is 
positive, which we assume to be the case. 

Q: It is unless you are on the top of the effective potential. 

A: Yes, but if you are here (a), you have to be careful about what you calculate. Usually, 
one wants to expand from here (b), i.e., around the local minimum. But I shall not 
discuss symmetry breaking in these lectures (except a little towards the end). 

Q: You mentioned that you need to put fi equal to T in perturbation theory. But I think 
in the end a physical quantity does not depend on /i... 

A: Yes, that is a very good point. Indeed, a physical quantity does not depend on /i. 
But remember that if you do the calculation in perturbation theory, the physical 
observables are independent of fi only up to terms that are of the same order of 
magnitude as those that are explicitly neglected. For example, in the case that I 
treated, I considered terms up to, and including, order A^, and the fi dependence 
will be of order A^. So at any finite order, there will be a residual /x- dependence. 
What I have argued is that one can exploit this dependence in order to improve the 
(apparent) convergence of the perturbative expansion. Such a strategy is sometimes 
referred to as the "principle of minimal sensitivity" , or "principle of fastest apparent 
convergence". At finite temperature, the coefficients of the expansion in powers of 
the coupling constant contain typically logarithms of the ratio fi/T. In order to avoid 
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these logarithms to become too large (and hence spoil the apparent convergence), it is 
judicious to choose /i ~ T (and detailed calculations suggest the more specific choice 
~ 27rT). 
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III. LECTURE 3 

Summary of last lecture 

Let me start by a brief summary of what we learned in the last lecture. 

What we did last time was to consider a scalar field with the hamiltonian density given 

by 

^ = + li^^^f + l^'^' + y{^) , (107) 

where we explicitly wrote V{(f)) = A0^/4!. Then we applied the general formalism to calculate 
thermodynamics, essentially the partition function, Z = Tre~'^^. Remember that this could 
be written as an integral over field configurations each configuration being weighted 

by the factor e"'^^'''^', with Se the Euclidean action corresponding to the hamiltonian H. 
The field configurations included in the integral are periodic in imaginary time. Here the 
calculation was organized as a power series, and we calculated the grand potential as 

n = log Z = Qo + X^i + A'/'fi3/2 + A'l^a + • • ■ ■ (108) 

In doing this calculation, we encountered two types of difficulties, namely divergences of 
some momentum integrals. Remember that loop integrals in Feynman diagrams involve the 
sum over the Matsubara frequencies Un and the integral d^k. 

We met two types of divergences. Ultraviolet divergence - these are familiar in field 
theory and can be handled by the general procedure called renormalization. In doing so, we 
introduced the notion of the running coupling constant, with the /3-function which describes 
how A varies with the renormalization scale n, 



Physical quantities, like the pressure, should be independent of /i, that is, 

f^^O. ,110) 

Of course, this holds if Q is calculated exactly. If you do an approximation on Q, for 
example if you calculate Q up to order A", Eq. (110) will hold up to corrections of order 
A"+^ which are not included in the calculation. In that case some residual dependence on 
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the renormalization scale will subsist in Q. This residual dependence can often be used as 
an indication of the accuracy of the calculation. 

Ultra-violet divergences are common in all field-theoretical calculations and they have 
nothing to do with the temperature. The formalism to handle these is well-established. 
We encountered also another type of divergences, called infrared divergences. These are 
intimately related to finite temperature effects. At the end of the last lecture I showed you 
that, in some particular cases, these infrared divergences are eliminated by performing an 
infinite resummation of Feynman diagrams, which corresponds to taking into account the 
generation of a "thermal mass". We've seen also that this resummation is responsible for 
the fractional power of A that occurs in weak coupling expansion of the pressure, namely, 
the term proportional to A^^^. Such a term is completely unexpected from the point of view 
of perturbation theory because if you expand naively the partition function in powers of 
A, you will only generate terms with integer powers of A, not A^/^. A term such as A'^/^ 
can only be obtained from resumming an infinite number of Feynman diagrams. Finally, 
we understood that there was one particular sector which is important, the sector with the 
Matsubara frequency Un = 0. Since Un = 27rnT, as soon as n is equal to 1 or bigger than 1, 
Un = 2TxnT provides an infrared cut-off in the propagator + fc^); therefore there is no 

infrared divergence for these non vanishing Matsubara frequencies. This ends the summary 
of the main issues that we discussed last time. 



My purpose today is to show you, from a very general perspective, why we have these 
specific difficulties at finite temprature and give you hints about how you can get around 
these difficulties by techniques that I shall explain in the next lecture. 

Before I do that, let me summarize the results that have been obtained for the pressure 
as a function of the coupling constant. To do so, let me change the notation slightly, and 
set = A/24 (then A'^^^ becomes ~ g^). The pressure P of 0^ theory is known up to order 
g^ (in fact, up to order g^hig [15]) 



The result is plotted in Fig. 14, divided by the free gas pressure Pq = (7r^/90)T^ (this has 



A. QCD perturbation theory and its breakdown 



P = Po + a2g^ + a^g^ + ■■■ + a^g^ + ■■■ 
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FIG. 14: Weak-coupling expansion for the pressure normalized to that of an ideal gas as a function 
of g(2'iTT) in (f)'^ theory. Taken from the review by Blaizot, lancu and Rebhan [9]. 

been calculated in the first lecture), and as a function of the running coupling constant 
defined at the scale 2ttT. As I argued last time, this is a natural scale to choose in the 
calculation of the thermodynamical functions. 

As the plot clearly indicates, perturbation theory does not appear to be very predictive 
unless the coupling constant is very small, g < 0.3. This same feature is also met in more 
complicated theories at finite temperature, such as QCD. In such theories, we know that 
the coupling constant decreases with the temperature, but unless the temperature is very 
large, so that the coupling is very small, the same pattern of bad apparent convergence is 
revealed. These calculations of high order contributions are technically demanding, but the 
result is disappointing. One of the questions that we have to answer is why it is so bad. 

1. Breakdown of perturbation theory 

In fact, before I turn to general considerations which will shed light on this question, 
let me show you that in QCD, we can reach a point where perturbation theory completely 
breaks down, whatever the strength of the coupling. Let me show you one particular example 
of a class of Feynman diagrams where the problem in question manifests itself. 

I shall consider the particular family of diagrams displayed in Fig. 15. The wavy lines 
represent gluon propagators. With the experience that we gained last time, we can expect 
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FIG. 15: Gluon loop diagrams contributing to the QCD thermodynamic potential. 

infrared divergences to occur when a;„ = 0. Since I am interested in these divergences, I am 
going to assume that a;„ = in all the propagators. Let us then examine the contribution 
to the pressure, or to the thermodynamics potential, of the £-loop diagram. For each loop, 
there is a factor of T that normally accompanies the sum over the Matsubara frequencies. I 
shall keep only the term n = 0, so there is no sum, but the factor T remains. And there is 
an integral d^k/(2TT)^ over the momentum (the factor (27t)^ is not relevant here and I shall 
drop it). Now there are also vertices. It is not hard to count their number: If there is one 
loop, there is no vertex; if there are two loops, there are two vertices; each time I add one 
loop, I add two vertices. Therefore, the power of g is 2i — 2. Now, remember that in the first 
lecture I emphasized the fact that the three-gluon vertex carries a momentum. Therefore, 
at each vertex is also attached a momentum k. Then there are propagators, of the form 
l/(a;^ + /c^). I set Un = 0, but I shall add a term which is a ficticious gluon mass, so 
that the propagator is 1 / {k^ + m^) . The number of the propagators has to be counted. This 
is 3{i — 1). You can understand this number in the following way: Start with two loops 
(second diagram in Fig. 15), there are three propagators . Each time you add a loop, you 
add two propagators. I can therefore write, very schematically, and for £ > 1) 



I am interested in what happens when all momenta are going to zero at the same rate. In 
this situation, all individual momenta can be combined to form a big momentum, K, in a 
space of 3i dimensions. Then, I separate the integral into an angular integral and a radial 
integral, and focus on the radial integral. In other words, in the same way as d^k = k'^dkdQ^ 
with solid angle dQ^, I'm going to write the integral as f {d^kY — )■ J K^^^~^^dKdQ-se, and 
drop the angular integral because, from the angles, no trouble is going to come. Let me do 
that. Then I get 




(112) 




(113) 
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i < 3 well-behaved integral 

£ = A T'^g^ log(r/m) 

e > 4 /T^ {g^T/mY"^ 

TABLE I: Order of magnitude of the £-loop contributions to the pressure. All the £-loop diagrams 
with £ > A are of the same order of magnitude when m = 0{g'^T). 

This integral is infrared divergent when m — )■ 0. To see that more clearly, let me replace the 
denominator + m? by K^, and put m as a lower cutoff in the integral: 

- T'g'''-'' [ dK-^^ . (114) 

If ^ < 3, this integral is perfectly well-behaved because the denominator has a vanishing 
or negative power of K. No problem. If £ = 4, we have a problem - we have an integral 
dK/K. There is a logarithmic divergence and I have put a cutoff m. In order to fix the 
scale there has to be another factor: the natural scale here is T . Thus, if I was doing the 
calculation at four loop, I would expect the final result to be of the form 

figcS - TV log(m/T) . (115) 

You see this is a divergent result: If I let m go to zero, which I should do because the gluons 
don't have any mass, then I would get an infinity. 

Let's look at £ > 4. This is even worse because, for £ = 5 for instance, I get J dK/K"^, 
which is a power divergence ~ 1/m. More generally, for £ > 4, we have 

^QCB-T'g^^'-'^ I dK/K'~' ~ T'g^^'~'\l/mY-^ = g^T' {g^T /my-^ . (116) 

J m 

This situation is not quite the same, but is very reminiscent of the situation that we met last 
time. I showed then that there were Feynman diagrams which where infrared divergent, the 
higher the number of loops, the more severe the divergence. Now, in the case of the scalar 
field theory, you know that there is a cut off because, as I showed you last time, the modes 
of the scalar field acquire a thermal mass of order of gT. Remember m|) ~ AT^ ~ g^T"^ 
(with A ~ g"^). You see, if the mass is of the order of gT, there would be no problem because 
g'^T/m ~ (7 ^ 1 in the weak coupling regime. In other words, the thermal mass cures the 
potential infrared divergences, and makes the successive loop corrections proportional to g^. 
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In the case of QCD, not all the modes acquire a mass of order gT. Only the electric modes 
do (the thermal mass is then identical to the Debye mass related to electric screening). But 
there are also modes in QCD which are magnetic in nature, and these modes may develop 
a mass, so-called "magnetic mass" . However, this magnetic mass is expected to be of order 
g^T. Then, because g'^T/m = 0{1), all the terms in the expansion (116) end up being of 
the same order of magnitude! This is a situation where there is not much you can do with 
a standard weak coupling expansion: if the mass is of the order g'^T, then all the terms in 
the perturbative expansion are of the same order of the magnitude. 

In the rest of this lecture, I would like to give you a simple understanding for why such 
nasty things happen. 

Q: Is there any different kind of mass other than magnetic mass in QCD? 

A: There are two types of modes in QCD, related to oscillations of the electric field and the 
magnetic field, respectively. The electric field behaves like a scalar field and it acquires 
a mass which is related to the screening phenomenon about which I will say more in 
the next lecture. That screening mass is of oder gT. So, there is no problem with 
the electric modes. The difficulty in the QCD plasma comes from the long wavelength 
magnetic modes. Such modes also exist in electrodynamics (a static magnetic field is 
not screened). But in QCD the magnetic modes interact with themselves, and it is 
this interaction which it is hard to calculate. 

B. Expansion parameter at finite temperature 

i. Harmonic oscillator at finite temperature 

The above is the review of the difficulties. Now, I will try to give you an insight into 
the physical origin of the difficulty and go through an elementary discussion of thermal 
fluctuations in quantum mechanics, going back to a very simple system, namely the harmonic 
oscillator, which you all know for sure. Basically, what we are doing with field theory is 
playing with an infinite collection of harmonic oscillators. You will see, this is a detour, but 
a quite instructive one. 
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Let's consider a harmonic oscillator in one dimension whose hamiltonian reads 

H = -p' + -u^x^ (117) 

In any state at T = 0, the expectation value of is equal to that of cu^x^ and given by 

(p2^ = ~ huj . (118) 

I shall drop the factor h in what follows (i.e., I shall use natural units where ^ = 1). Thus, 

(j92)~u;, {x^)^l/u. (119) 

Now, I need to do a little bit more elaborate things. Remember that, for the harmonic 
oscillator, it is useful to introduce creation and annihilation operators such that 

[a,a)] = l. (120) 

In terms of a and a\ we can write 

X = —^=[a + a' j , 



2uj 

1 1 2 

x^ = — (a + a^)^ = — (a'^ + a^ +aa^ + a^a), (121) 
2u 2uj 

and the eigenstates of H are of the form |no) = (a^)"^" |0). When we calculate the expec- 
tation value of x'^ in a given eigenstate of the harmonic oscillator, e.g. (n| \n), we get 
contributions only from the last two terms. Using the commutation relation, one easily gets 

(nolxVo) = 7^(l + 2no). (122) 

Here you must recognize something: the factor (1 + 2nk)/{2uk) is indeed that same as that 
in the integral giving the fluctuation of the scalar field (see for instance Eq. (87)). 

Now, at finite temperature, T 7^ 0, we can apply the same results. On the average, states 
will be occupied with a probability given by the Boltzmann factor, and the expectation value 
of at a given temparature T will be given by a formula analogous to Eq. (122), namely 

{x')^r^^{l + 2nT), (123) 

where 
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There is a limit which is interesting to us. It is the high temperature hmit, T ^ co. 
In that case, u/T is a small number, and we can expand the exponential in the statistical 
factor ut- In this regime, Ut ^ T/u ^ 1, and we can ignore the contribution of the zero 
point fluctuations (the "vacuum part"). What remains is 

T 

u 



^")t - -2 ■ (125) 



This result is familiar and reflects the equi-partition of the energy, 

as expected from Boltzmann statistics. 

Of course, the harmonic oscillator itself is fine, but what I'm really interested in is to 
study the effect of the interactions. So let me add interactions by changing the hamiltonian 
H ^ Hq + Hi. What's a good choice for iJi? To stay as close as possible to the 0^ scalar 
field theory, I take 

Hi = ^x' . (126) 

We may now calculate the energy levels as a function of A which is supposed to be a small 
parameter. So, the energy levels will be written in the form 

Eno = {no + ^)uj + aiX + a^X^ + ■■■ . (127) 

where the numbers Oj depend on uq. At finite temperature, the pressure, or the thermo- 
dynamical potential Q, are functions of A, and we are interested in the expansion of these 
functions in powers of A. Then, the question we want to ask is what controls the expan- 
sion. We want Hi to be small compared to Hq in some sense. Let this be measured by a 
parameter, 7: 

7 - ^- (128) 

At this moment this is not very well-defined. Let me specify this parameter a little bit better 
by taking some average values: 
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where I assume {x^) ~ (x^ y. This 7 is my dimensionless expansion parameter. As you see, 
this depends on the strength of the couphng A, and obviously A should be small for the 
expansion to make sense. But it also depends on (x^) and w^. 

At T = 0, what is 7? At T = we have calculated (x^) ~ l/w, so that 

7-4- = 4- (130) 

You can verify that the dimension of A is the same as that of w^, so that A/w^ is dimensionless, 
as it should. If you are doing a perturbative calculation of ground state properties at zero 
temperature the calculation is valid provided that \/uj^ is small compared to one. 

At T 7^ 0, things are different. In the high temperature regime, T ^ cj, we have 
(x2)^ ~ T/u;2 and 

i--A-')--A = \-- (131) 

The factor X/u^ is the familiar one, but there is another factor which can be big when the 
temperature is big. So, at finite temperature, it's not enough to have X/u^ ^ 1. If the 
temperature is big enough, 7 may become of order unity, even if X/u^ <^ 1. 

What happens physically at finite temperature is that the system expands, because states 
with large quantum numbers become occupied, leading to an increase of (x^). The situation 
here is comparable to that of perturbation theory for the excited states at T = 0. There 
the expansion parameter would be 7 ~ {X/uj^)no, which may become of order unity if no is 
large enough. 

Remark. One way to improve perturbation theory is to define an effective frequency as 



^ujj^ = uj' + i^X{x') , (132) 



that is, absorb part of the interaction into the "unperturbed hamiltonian" , i.e., replace 
x^ — )■ x^(x^), and adjust (x^) self-consistently. This leads to an approximation which is 
similar to the Hartree approximation of many-body physics. 



2. Field theory 

Let's move now to field theory. The field theory is almost identical except for one impor- 
tant feature, namely that a mode carries a given wavelength. We are thus led to distinguish 
between layers of fluctuations at different wavelengths. 
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Q: I want to interrupt to make sure that this difficulty comes from some kind of divergence 
at small u. 

A: In the harmonic oscillator there is no divergence because we have only one mode. The 
difficulty comes from the fact that, if the temperature is big, i.e., T ^ w, a large 
number of quanta are excited, and this may make the expansion parameter of order 
unity even if the zero temperature expansion parameter A/w^ is small. If X/uj'^ = 10~^ 
and ii T/u = 100, you are expanding in powers of lO's — that is the problem. 

Q: This factor comes from the singularity of the distribution function? 

A: It comes indeed from expanding the distribution l/(e"/^ ~ !)• In the regime where the 
temperature is huge compared to the distance between the energy levels, we need to 
take into account many levels, all of which contribute to produce a large value of (x^). 
(A similar difficulty would occur if you were calculating corrections to the energy of 
an excited state with large uq.) 

Q: At high temperature the system becomes classical. Is that the reason why you can use 
this mean-field approximation to treat the higher-order terms? 

A: Well, this Hartree approximation that you are referring to can be used also at zero 
temperature. It turns out that this is a rather good approximation, even though it's 
more useful at finite temperature. But I don't think that this is related to the classical 
approximation. Of course, if you are in the classical regime, then (x^) is big, and it is 
advantageous to do this. 

Q: I thought that the quantum fluctuations become less and less important ... 

A: Yes, that's right. In the expectation value, (x^) ~ 1 + 2nT, the ffist term, which is 
the vacuum fluctuation can be ignored if ut is large enough, which is the case at high 
temperature where ut ^ T/u ^ 1. 

So let's move now to the field theory. We replace x"^ by 0^, and (0^) is now given by an 
integral 




(133) 
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where Uk = l/(e'^'=/^ — 1). I am also going to assume that m = so that = k. In the last 
(approximate) equality, I have just ignored the vacuum fluctuations. Compare this to the 
formula (122) that we had before for the harmonic oscillator. The main difference is that 
we have now a sum all over the modes which are labeled by the momentum k. 

Now, I want to repeat the analysis of the expansion parameter; I want to compare "the 
kinetic energy" and the "the potential energy." The kinetic energy is 



;(V0)^> , (134) 



and the potential energy is 



A (04)^/ (02)^ , (135) 

where I assumed that {(p^) ~ (0^)^, which is good enough for the qualitative discussion that 
I want to present. 

Now, comes an important remark. The integral in Eq. (133) is dominated by the large 
momenta. One way to see that is that, in the absence of the statistical factor, it would be 
quadratically divergent. What the statistical factor does is to provide a cutoff at k = 0(T). 
Let us then consider the integral with an upper cutoff k: 

. r#^^. (136) 

J (27r)3 cfc ^ ^ 

This integral, whatever k, is, is always dominated by the largest possible momenta, and 
therefore what will contribute dominantly in this integral are the fluctuations which have 
momenta k of the order of k. I'm going to call (0^)^ the "contribution of fluctuations at 
scale k" . In the same spirit, I'm going to define the typical kinetic energy in (134), at the 
same scale k, as 

((V0)^>^~K^(0^>^. (137) 

Now, I have all the tools to define the expansion parameter 7, as we did in the analysis 
of the harmonic oscillator. The big difference is that now this depends on n. I call 7^ the 
ratio the potential energy ~ g'^ (0^)^ to the kinetic energy (0^)^) that is, 

7. ~ - m . (138) 

The rest of the discussion will be concerned with the analysis of the typical momentum 
scales that appear in the ultra-relativistic plasmas in the weak coupling regime. 
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hard soft ultra-soft hard-soft* 
g'^T gT-T 

1 1 

TABLE II: Expansion parameter for modes with typical scale n. *see text for definitions. 



K T gT 

Ik g 



C. Interplay of the various wavelengths 

1. Self-coupling of hard modes 

There are "natural" scales in the ultra-relativistic plasma; one scale, which I call "hard" , 
is K of the order of T. This is the scale for the typical plasma particles which are the modes 
of the field that carry momenta of the order of the temperature T. Remember that the 
density of such particles goes like T^, and the average distance between them is therefore 
of order 1/T. Note that 1/T is also the typical de Broglie wavelength, so that these hard 
particles are always quantum ( the photons in the blackbody radiation). 

Then I want to calculate (0^)^, with k ~ T. This is an integral that we have already 
done. We get 

7T ~ ~ . (139) 

For these particular modes, there is therefore no problem: The expansion parameter is g^, 
so that perturbation theory is OK, provided that g is small. (In fact, if you refer to a 
familiar quantum field theory like QED, you know that the expansion parameter is g'^/'iir 
rather than g"^. ) 



2. Self-coupling of soft modes 

In an ultrarelativistic plasma, there is another scale, which is called the soft scale, k, ~ gT. 
You will see why this particular scale emerges dynamically. We already have an example of 
that scale with the thermal mass, which is proportional to gT. 

When K = gT ^T, 
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so that 




9T (^r)2 




(141) 



You see an interesting phenomenon here. You see that 7^ is still a good expansion parameter, 
because it proportional to the coupling constant, which is a small number. But you also see 



the self-interaction of the soft modes is still valid, but the expansion is not in powers of 
but in powers of and so is less precise. This new expansion parameter, g rather than g"^, 
where g is the strength of the self-interaction between the soft modes, is the origin of the 
odd powers of g that we have encountered in the expansion of the pressure. 



But there is something more, which is very important. The parameter ~ g compares 
the strength of the self interactions of the soft modes with their kinetic energy. Now, you 
may ask another question which is as relevant as this one. You may ask what about the 
possible coupling between the soft fluctuations and the hard fluctuations. 

Kinetic energies of the soft modes are of order (gT)"^, while the interaction energy between 
the soft and hard modes is ~ g^ {<P'^)t ~ 5'^^^- So the motion of a soft mode, with wavelength 
of order l/{gT), is non-perburbatively renormalized by its coupling to the hard degrees of 
freedom. Here you cannot expand because the two effects are comparable. So, this is a non- 
perturbative correction, which goes under the name of Hard Thermal Loops. The thermal 
mass which I have introduced the last time, and also discussed today, is one particular 
example of this correction. 

So you see, at the soft scale the landscape complicates a bit. We know that we can treat 
the soft modes perturbatively, although the perturbation theory is not as accurate as for the 
hard modes. But the propagation of the soft modes itself is affected by their coupling to the 
hard modes in a non-perturbative fashion. In order to treat this phenomenon, perturbation 
theory is not enough. But this is well under control. The mass resummation that we have 
discussed last time is a simple example of what needs to be done to handle this problem. 

Q: I thought that you defined the 7^ as the ratio of the kinetic energy to the potential 
energy of each mode at the scale k, right? 



that, while we had 77- = g 



^ before, now 7^7- = g. It means that the perturbation theory for 



3. Coupling between the soft and hard modes 
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FIG. 16: Schematic representation of the expansion parameter for the soft modes k = gT and 
the couphng ^qT-t of the soft modes to the hard modes. 

A: Yes, this is what I've done before. But here I'm defining another 7. I look at a mode 
with wavelength 1/gT, much bigger than the wavelength 1/T of the typical plasma 
particles. When such a mode propagates, it can interact with the hard modes through 
loop corrections (see Fig. 16). The present expansion parameter is coming from the 
comparison between the kinetic energy of the soft modes, which is {gTy, and the 
potential energy coming from the hard loop, ~ T^, times the coupling g"^. I'm arguing 
that they are of the same order of magnitude. 



4. Ultra-soft modes 

Now let me move on to the another natural scale, and confront a real catastrophy. This is 
the ultra-soft scale, k, ~ g^T. The reason why this scale occurs is that it is at this particular 
momentum that you have the complete matching between kinetic and potential energies. 

At the scale k ~ g'^T, we have indeed 

(0^)^,^ ~ T X = g'T^ (142) 

so that 

(W)^>,.^ - {9'Tr{<P')^.^-9'T\ (143) 

Then you see that the kinetic energy and the potential energy are of the same order: 

Ig^T 1 • (144) 

No expansion is possible! In other words, the ultra-soft modes remain strongly coupled for 
arbitrary small coupling. Even if the coupling is 10~^^, you cannot expand in the powers 
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of 10 because the kinetic energy and the potential energy are always of the same order. 
Note that the contribution of the ultra-soft modes to the thermodjTiamic potential is 

g' {<P%r - 9\g'T'f = g'T' . (145) 

This is the order g^T^ at which QCD perturbation theory breaks down, as we have discussed 
earlier. 

You see that the difficulties with QCD perturbation theory at finite temperature can be 
understood from rather general considerations. The key point here is that ultra-relativistic 
plasma should be viewed as multi-scale systems. These are difficult to treat because, typi- 
cally, approximations devised for one particular scale do not work for other scales. 

Perhaps you have heard about "the strongly coupled quark-gluon plasma" and about 
the AdS/CFT correspondence which allows us to calculate at infinite coupling. From the 
perspective that I have just outlined, the problem with finite temperature filed theory is not 
so much connected with the absolute strength of the coupling. The picture provided by the 
AdS/CFT correspondence is one in which all modes interact with infinite strength. What 
I have tried to argue is that this picture ignores the important fact that the effects of the 
interaction in the QCD plasma depend on the wavelength of the modes that one considers. 

Q: This discussion reminds me very much of the critical phenomena. Near the critical 
phenomena, we have long wavelength massless modes. They interacts very strongly 
and there are many scales intermingled. Is there any insight to be gained from the 
experience of the critical phenomena? 

A: Well, in fact, in the last two lectures, we will be dealing with one particular critical phe- 
nomenon, which is Bose-Einstein condensation. This very same problem will appear. 
The techniques to solve it will be borrowed from what I will shall discuss next time. 
Here, as I said before, the ultra-soft scale does not occur for the ordinary scalar field, 
because a thermal mass is generated. But if you are in the vicinity of a second order 
phase transition, you can adjust a parameter in such a way that the effective mass at 
some temperature vanishes. Then we are exactly in the same situation. 

Q: I think the origin of the classical, long wavelength modes is very similar. 

A: You will see that the structure of the infrared divergences in the analysis that I will do 
quickly for the Bose-Einstein condensation is indeed very similar. That is why I put 
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the two topics together in these lectures. In the last lecture, I hope to have some time 
to introduce you to some of the modern techniques of the renormalization group to 
handle explicitly such multi-scale phenomena. 

Q: Can the breakdown of perturbation theory at order be overcome by using other 
weak coupling techniques? For example, the mass-screened perturbation theory, hard 
thermal loop resummation, optimized perturbation theory, or others? 

A: I have discussed here the origin of the difficulty. I have not discussed at all what we 
should do in order to overcome it. We know that if you use perturbation theory in a 
regime where soft momenta are integrated over you should not be surprised to meet 
infrared divergences because the expansion that you are trying to use makes no sense. 
That we know a 'priori without doing any calculation. So we know a priori that we 
should do better. 

Now what to do depends very much on what you want to calculate. I will briefly 
comment next time on the application to the thermodynamics of QCD. For the ther- 
modynamics, hard degrees of freedom are dominant. So you expect that perturbation 
theory, plus some corrections, will work, and indeed it does. So, you can use screened 
perturbation theory, you can use the 2PI formalism, you can do hard thermal loop 
resummation, all these techniques will essentially deal with this part of the problem. 

If you want to address the situation described here, and if you want to calculate ex- 
plicitly the contribution of the very long wavelength modes (the modes that I called 
ultra-soft), then there is no other way, that I know of today, than to do lattice calcu- 
lations. 

But for the thermodynamics, observe that the long wavelength, small momentum, 
modes have very small phase space, so their contribution to the energy density and 
the pressure is presumably small. But if you are asking about correlations at long 
distances in the plasma, or perhaps transport phenomena, you may need to worry 
about these modes. Their analytical treatment remains an open issue. 

Q: You said that the coupling between the soft and hard modes becomes in a sense dan- 
gerous because the ratio becomes essentially one. If you consider the coupling of the 
ultra-soft modes to the other modes, what happens? 
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A: This is an interesting question. I let you meditate about it. I shall just point out that 
the scale gT is uniquely determined by the requirement that ~ 9'^{4''^) k{<P'^)t- 

Next time I will show you the technique of the effective field theory. This will be directly 
relevant to what I will do for the Bose-Einstein condensation. And I will also show you how 
the hard thermal loop emerges in a dynamical context, i.e., from kinetic theory. This is an 
interesting perspective, which may be also of special interest to those of you who are are 
working on kinetic equations. 
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IV. LECTURE 4 

Summary of lecture 3 

One of the things that I have emphasized in the previous lecture is that the quark-gluon 
plasma, or more generally an ultra-relativistic plasma, can be viewed as a multi-scale system: 
there is only one scale to start with, namely the temperature, but at weak coupling, other 
scales are generated dynamically. Such a plasma contains modes with various wavelengths, 
in fact there is a continuum of wavelengths. And the important point is that the effect of 
the interaction between these modes depends very much on their wavelengths. In order to 
characterize this interaction I have introduced the quantity (0^)^ which I referred to as the 
fluctuations at scale k: 

/ ,2\ d^/i; rife 1 

(Remember that the dominant contribution to this integral is determined by momenta which 
are of the order of the upper momentum cutoff n.) This allowed me to define an expansion 
parameter which controls perturbation theory, that I called 7^: 

Ik = 2 ~ • ^^^') 

The discussion of the last lecture can be summarized in the table below. We have consid- 
ered three particular scales, T, gT, and g^T. I'm assuming here that g is a small number, 
g 1. (This is a condition which in fact does not need to be that strict because there are 
factors 27r, for instance, that have been left out. ) 







g'{4>'). 


Ik 


K = T 




g'T' 


g' 


K = gT 




g'T' 


g 


K = g'T 


g'T' 


g'T' 


1 



The fluctuations at scale T have interactions controlled by g'^. For these, the perturbative 
expansion behaves as at zero temperature (and is really an expansion in powers of g'^/Air). 
If we move down to the scale gT then we have again an expansion parameter which is small. 
But it is of order g instead of order g"^. That means that perturbation theory for these 
particular set of wavelengths will be less precise. The expansion will not be in powers of 
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g'^ but will contain odd powers of g. The contributions of order g^ will be at the center of 
our discussion today. And finally at the scale g^T the expansion parameter is of order one, 
which means that the very long wavelength modes will remain strongly coupled, however 
small g may be. 

I also made another remark concerning the scale gT, which will introduce today's discus- 
sion. The remark is that if I compare the kinetic energy of a mode with momentum k ~ gT, 
with the contribution of its interaction with the fluctuations at scale T, I find that they are 
of the same order of magnitude, namely: 



This indicates that even though the self-interactions of the modes with momentum gT 
are of order g, the motion of these modes is strongly modified by their coupling with the 
fluctuations at the scale T. This feature is at the heart of what is known in the literature as 
the "hard thermal loops". And I will tell you a good deal about these hard thermal loops 
today. 

What I would like to do is to explain how we can handle this coupling between soft and 
hard degrees of freedom. This will be done by introducing an important construct which is 
that of effective field theories. I shall do that both in Euclidean, or imaginary time formalism, 
and also in real time, where the effective theory takes the form of a kinetic theory. This 
program would require a whole set of lectures in itself, which of course I cannot do in one 
afternoon, so I shall have to skip some details, especially in the second part of the lecture. 



Let me now discuss effective theory. What we have to deal with is a situation where 
degrees of freedom with different wavelengths are coupled together and interact differently 
depending on their momenta. I'm going to approach the problem within the path integral 
formalism. Let me remind you the formula we have for the partition function 




(148) 



A. Effective theory 




(149) 



where 




(150) 



59 



What I am going to do is a more elaborate version of something that I have already intro- 
duced. Remember two lectures ago. We discussed the situation where (3 goes to 0. If /3 goes 
to 0, and if nothing singular happens then, we can ignore the time dependence of the field. 
Then Se simplifies since the integration over time gives just a factor /3: 

S^^pJ d'x (1(V0)^ + + ^04^ . (151) 

Since P is the inverse of the temperature, this corresponds to the high temperature limit. 
In this limit, the quantum filed theory reduces to a classical three dimensional field theory. 

What I want to do now is essentially an elaboration of this remark. I will often refer to 
this high temperature limit as the "classical field approximation". As I have also argued 
earlier, this is the approximation that is obtained by ignoring the non vanishing Matsubara 
frequencies when one expands the field in Fourier space. That is, this is the approximation 
where one keeps only the vanishing Matsubara frequency component of the field. 

So let me consider the Fourier expansion of the field: 

= T0o(x) +tJ2 e~^"""</'n(x). (152) 

What is done in Eq. (151) is just keeping 0o and ignoring all the rest. What I would like to 
do now is to show you how we can take into account the rest. I do this in the framework of 
a general (p^ scalar field theory. (The procedure generalizes to more complicated theories, in 
particular to QCD, but we shall not have time to discuss this in detail.) This approximation 
will play a crucial role in the solution of the problem that I shall discuss in the next two 
lectures, namely Bose-Einstein condensation. 

How are we going to take into account the effects of the components with non vanishing 
cjn? Note that I can make a change of variables, going from a path integral over field 
configurations in spatial coordinates x time, towards configurations in spacial coordinates 
X frequency. In other words I can write the measure of the path integral as a product 
X Yln^o ^'^n- (■'■ concerned here just with the time dependence, and I do not write 
explicitly the dependence of the field on the spatial coordinates.) Then I can rewrite the 
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path integral as follows 



1.^ 



oe~'^°ff[*°], (153) 



where 



g — S'efl[<}!'0] 



1 JjD0„e-^^['^O''^"l . (154) 

How can we calculate S'cfr[0o]? In fact, you know enough to be able to calculate S'efr from 
what I have told you already. Look indeed at Eq. (149), and compare it with Eq. (153). You 
see that exp (—S'efr) is a partition function, for a system in which 0o appears as a frozen, 
given parameter. So we know how to calculate that. We know that the thermodynamic 
potential f2 in Eq. (149) is given by the set of all connected Feynman diagrams. Similarly, 

S'efr [0o] = set of all connected Feynman diagram with external lines "0o" (155) 

We know the leading order contribution to Se//[0o]: this is the "tree- level" action, ob- 
tained by evaluating the original classical action with the field 0(r, x) = T0o(x), and it is 
identical to Eq. (151). In fact it is convenient to rescale 0o by a factor -\/T, that is to define 
(pQ as 

(/>o(x) = / dr(/.(r,x) (156) 

With this definition, you see that a term such as /3</)^ in Eq. (151) would lead a term Tcf)^ 
with the definition of Eq. (152). The new definition absorbs the factor T into leaving 
the simpler expression for the leading order effective action: 

1 , ^2 ,2 AT 

2(V0o)' + ^0o + ^ 

Corrections to the leading order are generated, for instance, from a loop expansion in the 
theory for the hard modes, whose Euclidean action is given by 



Sfh\M = ld'x{ ^(V0o)^ + + (157) 



Se[(Po, ^n] = Y.) ^V0„V0_„ + ^ (m^ + XT^l + col) 4>n4>~ 



nmk nmkl 



where the fields 0„ have been rescaled by the same factor vT as 0o- For instance a one loop 
contribution is easily obtained from the part of SE[4>o,4'n\ that is written in the first line 
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FIG. 17: One loop dia gram with external lines (^O; and internal lines (pn- 

of Eq. (158), and which is quadratic in The one loop diagrams produce a correction to 
the mass (to be discussed later), a correction to the four-point function (see the Fig. 17), a 
contribution to the six point function, etc. There will be other corrections generated by the 
terms in the second line of Eq. (152), but these will be of higher order in A. For instance 
the first correction to the mass generated by the term with three 0„ will be of order at 
least. 

At this point I need to specify more carefully the separation between what is meant by 
"soft" and "hard". Naively we may attempt to call soft the mode with n = 0, and hard the 
modes with n ^ 0. This is what we have done so far, but this is not enough. If we just leave 
things at this stage, we are going to generate an effective action which will contain arbitrary 
powers of 00, but which will be mostly non local. For instance, consider the diagram in 
Fig. 17. As it stands, it is a complicated function of the momenta carried by the various 
lines labelled (pQ. We would like however to continue working with a local effective action, 
that is, we would like to be able to consider this diagram as a correction to the coupling 
constant. This can be achieved by introducing a separating scale A, chosen so that (here I 
rename X ^ g"^) 

< A < 2nT, (159) 

because I'm now assuming "soft" will involve some energy scale of order gT, while "hard" 
will involve some energy scale of order 27rT. Then I redefine the separation between soft 
and hard in the following way. In the soft sector I consider the mode n = 0, but I also 
assume that the momenta are limited by A. So, when I calculate an integral in the effective 
theory, I do all the momentum integrations up to the scale A. On the other hand, the hard 
sector contains all the modes with non vanishing Matsubara frequencies, plus the sector 
with the large momentum components {k > A) of the mode n = 0. We have now a cleaner 
separation. In all cases you see that, with = u"^ + k^, the separation which I have 
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introduced guarantees that the hard momenta satisfy K > K and the soft momenta satisfy 
K < A. 



hard 


n 7^ + n = and 


k 


>A 


soft 


n = 0,|fc| <A 



Thus defined, the effective theory for soft external momenta is vahd only when the momen- 
tum p carried by the field (^o is smaller than A. On the other hand, the loop integrals are 
dominated by hard momenta. 

Now, we have gained something with respect to the argument of locality that I alluded 
to earlier. We know that the momenta which are inside the loop are all going to be large 
compared to the momenta outside. The locality will result from the fact that one can expand 
mp/K where i^' is a typical momentum inside the loop, and p an external momentum. The 
expansion in powers of p/K means the expansion in the field and its derivatives. Therefore 
that means that the effective action is going to be given by a series of monomials built from 
the field 0o and its derivatives: 

C = a(j)l + b (V0o)' + + d(j)l + e (0oV0o)' (160) 

Note that this expansion preserves the symmetry of the original lagrangian under (po — )■ —(po. 
In principle we have an infinite collection of terms, and we need a guiding principle to 
truncate this expansion, otherwise we cannot do calculations. We shall use here a weak 
coupling approximation, where the successive terms in the effective action can be calculated 
using perturbation theory. 

Let me give you the form of the effective action with more standard notation: 

^eE = /a + / ( V0o)' + ImI<PI + |<^o + ^00 + ■ ■ ■ ) , (161) 

where /a represents the contribution of the hard mode to the thermodynamical potential. 
Note that the dimensions of the various terms are characteristic of a three dimensional field 
theory: 0o have mass dimension 1/2, [0o] = M^, has a dimension of mass, [g'^] = M, and 
h is dimensionless, [h] = M°. 

The partition function of the system can be written as 




(162) 
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If Seff is calculated exactly, this is identical to Eq. (149). The strategy now is to obtain Z 
from an approximate iSefrf^o]- The coefficients in Seff can be calculated using perturbation 
theory. Some are easy to get, because they exist already at tree level. This is the case for 
instance of gl (the corresponding tree level diagram is given by Fig. 18, whose value can be 
red off Eq. (157): 

gl = g'T (163) 

The next vertex, h, has no contribution at the tree level because there is no term like (p^ in 
the original lagrangian. The term of order 6 in 0o is induced, at leading order, by the one 
loop diagram of Fig. 19. This diagram is proportional to g^, and is therefore subleading if 
g is a. small number. 

Everything is fine so far with the general strategy. We have a very systematic way of 
calculating the contributions to the effective action. One can do that in perturbation theory, 
i.e., one can calculate the coefficients of the effective action from Feynman diagrams. Note 



64 




FIG. 20: Tadpole diagram that contributes to the coefficient of (/)q in the effective theory 

that there is no contradiction here with the fact that perturbation theory cannot be used to 
calculate within the effective theory. I am using perturbation theory here in the sector where 
it is safe, because in calculating the coefficients of the effective theory, I am only integrating 
over the hard modes. To calculate with the effective theory, I have to do something more 
sophisticated. 

B. Calculation of the thermal mass 

What I shall do now is perform a simple calculation using the effective theory, in order to 
illustrate how things work. I shall calculate the thermal mass. The calculation will proceed 
in two steps. First, I shall calculate explicitly the one- loop correction to the coefficient of (^q 
in the effective action (161). Then, I shall use the effective theory to calculate the correction 
due to the soft modes. One issue that I want to address is that of the dependence of the 
results on the arbitrary scale A which comes in as soon as we consider loop corrections. I 
shall verify that when the calculation of a physical observable is correctly performed, this 
dependence disappears from the final result. 

1. Contribution of the hard modes 

The coefficient of 0q in the effective action is given by all the Feynman diagrams which 
have two external 0o lines. At leading order, there is a single diagram, often called the 
"tadpole" diagram, displayed in Fig. 20. This is easy to calculate. We have already done 
so a number of times. But now, we have to pay attention to the fact that the loop integral 
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runs over hard momenta only. Thus, we have 

M^w f d^k 1 - 6no + Snoejk - A) 

M(A) = -T>_^y^ ^''^^ 

where I have taken into account that the hard degrees of freedom are all the modes with 
n ^ (hence the contribution proportional to 1 — 5„o) ? as well as the contribution from the 
mode n = with momenta bigger than the dividing scale A (hence the term 6no0{k — A)). 
Now, use the following relation between 6 functions: 

^(A; - A) + ^(A - A;) = 1, 
l-Sno + Sno (1 - ^(A - A;)) = 1 - 6noe{A - k), 

and rewrite the integral in (164) as 

.2 r ^3/ 



2 (27r)3a;2 + p 3 J {2tiY k^' 

The interpretation of the first line of this equation is simple: the term 5nod{A — k) explicitly 

removes from the loop integral the contribution of the soft momenta, which needs to be 

calculated more accurately. This will be done by using the effective theory. By performing 

the sum over the Matsubara frequencies, one obtains 

,.2/»x g'^rr. f d^k l + 2nk d^k 1 , , 

A = / - — - , - — T / - — - — . 166 
^ ^ 2 J (27r)3 2A; 2 J (27r)3 P ^ ' 

As we have done earlier, we shall drop the term which is ultraviolet divergent and indepen- 
dent of the temperature (this goes away with a standard ultraviolet renormalization of the 
mass at zero temperature). The result can be then written as follows 

9^ / _ A_6 
24 V Tn' 

Since A is much smaller than T, the second term of (167) can be viewed as a correction, but 
which depends on A. As we shall see, this A dependence will cancel against an analogous 
contribution from the effective theory calculation. 



2. Contribution of the soft modes 

What I have done so far is to calculate the coefficient in the effective action. What 
I want to do next is to calculate the correction to the thermal mass in the effective theory. 



^'(A) = ^(l-;^-^), A«r- (167) 
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FIG. 21: Tadpole diagram in the effective theory. The loop involve the propagator of the field (j)o 

What is the correction to the mass? This is the correction that results from the self interac- 
tion of the field 0o- The correction is given by the same diagram as in Fig. 20, but redrawn 
in Fig. 21 in order to emphasize the elements of the calculation. The diagram in Fig. 21 is 
a diagram in the effective theory. So the vertex is g^. And the loop integral involves soft 
modes only, that is, the momenta are limited to A. As for the propagator, it's inverse can 
be read off the effective action itself in Eq. (161): k"^ + M^. The calculation then proceeds 
as follows 



1^ lo k?+Ml 



(27r)3 fc^+Mj 



1 _ Ma tan-1 

^ A Ma 



:i68) 



Remember that Ma is of order gT (see Eq. (167)), so that A ^ M\ ~ gT. By expanding 
the second term of (168) for large A/Ma, and keeping only the leading order correction in 
Ma, that is, M| ^ M| = g^T'^/24: {gl = g'^T), one obtains 

~ - ^MT. (169) 

It follows that the physical mass is given by 



24 V Tn^J 47r2 Svr' 
g'T' g'T gT 



(170) 



24 Stt 

You see that, as anticipated, the A dependent term coming from the parameter of the 
effective theory. Ma, and that coming from the cutoff in the loop integral within the effective 
theory, precisely cancel. The correction, of order g^T^, agrees with the calculation done in 
the second lecture (see Eq. (105)). 
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I have presented a very simple example of the construction of the effective theory at the 
scale gT. In the particular context of QCD, such an effective theory (also called dimensional 
reduction), has been pushed to a high degree of accuracy. When I mentioned earlier in the 
lecture the calculation up to or g^logg for the QCD pressure, these were obtained by 
relying on such techniques (see for instance [16] and references therein). Of course in QCD, 
these techniques are more elaborate. In a gauge theory, you cannot simply put a cut off on 
integrals, as we did. You have to use more sophisticated regulators. The technicalities are 
more difficult to master. But the basic concepts can be understood from the simple scalar 
field theory discussed in this lecture. 
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C. Real time Hard Thermal Loops 

In this second part of the lecture, I would like to give you another perspective on hard 
thermal loops, using kinetic theory. As we shall see, kinetic theory emerges as the effective 
theory that allows us to efficiently handle the coupling between hard and soft degrees of free- 
dom in ultrarelativistic plasmas. I shall not proceed through a systematic derivation, which 
would require more lectures, but shall try to indicate the main steps in such a derivation 
and emphasize the main physical aspects. 

1. Real time propagators 

As a preliminary, I would like to comment about connections between the imaginary time 
and the real time formalisms, and in particular remind you of some relations based on the 
analyticity of the propagators. In real time, we define 

G{t^M) = {mtMt2)) = |Tr[e-^^r(0(ti)0(t2))] (171) 

(T0(ti)0(t2)) = e{ti - t2)(0(tl)0(t2)) + e{t2 - tl)(0(t2)0(tl)) (172) 

where I have omited the spatial coordinates (which play no role in the discussion) in order to 
alleviate the notation. We also define G^(ti, ^2) = (0(^1)0(^2)) andG<(ti,t2) = (0(^2)0(^1))- 

Let us focus on G^(ti,t2). By making explicit the time dependence, and expanding on a 
complete set of eigenstates of the hamiltonian, we get 

G>{ti,t2) = -Tr(e-^^e'^*i0e-*^*ie'^*20e-*^*2) 

= ^^e-^^"e*(*l-*2)^"e-*(*^-*2)^™|(r^|0|m)|^ (173) 

n,m 

or, setting ti — t2 = t, 

G>{t) = ^^{e-''^"e'^"V^'"*|(n|0|m)|'}. (174) 

n,m 

In this expression, we can, as we have already done several times in another context, set 
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FIG. 22: Analyticity domain of G>{t) 

it = T oi t = —IT. We then get 

G>{-it) = 45"{e"^''"e"^"e-"^™|(n|0|m)|'}. (175) 

n,m 

Let us see under which conditions this substitution is legitimate, or more generally, under 
which conditions the time can be given an imaginary part. In most cases, the convergence 
of the sum will be controlled by the exponential factors. If r > /3, e'^^" will win compared 
to e~^^" and the sum will explodes. However if r < /3, e~^^^ dominates and this will kill 
the other terms as En gets large. Therefore the sum over states in Eq. (175) that allows 
the calculation of G^{t = —it) is finite if r < /3. One concludes that G^{t) is an analytic 
function of t in the strip —i[5< \m.t < 0, as indicated in Fig. 22. 

After Fourier transform, the propagator enjoys also analyticity properties in the frequency 
complex plane. Recall the form of the propagators in terms of Matsubara frequencies 

uj^ + k'^ + rn? —{iUnY + k'^ + m'^ ^ ^ 

By changing iujn to I transform this into the familiar propagator of a relativistic particle 

1 1 



— . (177) 

— {iujn)'^ + /c^ + m? —uj"^ + k"^ + rn? 

Once you have continued the propagator from discrete imaginary frequency to an arbitrary 
complex frequency uj (see Fig. 23), you make apparent the pole at a real frequency Uk 
corresponding to the excitation energy of the system with plus or minus one particle, Uk = 
±v^A?"+~?T?. Thus, the analytic continuation of the Matsubara propagator allows us to get 
information on the excitation energies of the system. 
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FIG. 23: Analytic continuation in the complex lo Plane 

However, it is sometimes difficult to do the analytic continuation exphcitly. For instance 
you may know the propagator only numerically, at all values, or only at a given subset of 
values, of the Matsubara frequencies. Then, performing the analytic continuation in order 
to extract the physical singularities may be an (almost) impossible task. Therefore, it may 
be advantageous to be able to perform calculations directly in real time. This is what we 
shall do in this lecture. 

2. An example of hard thermal loop 

After this reminder about analyticity property, I would like to show you one particular 
Feynman diagram calculation of a hard thermal loop. This is actually how hard thermal 
loops were discovered, by explicitly calculating a series of Feynman diagrams, and making 
the appropriate kinematical simplifications. I shall do a calculation of a one-loop self-energy 
in a scalar theory with a 0^ interaction, in order to be able to compare with corresponding 
results in electrodynamics that I shall consider next. The scalar 0^ field theory is not 
completely stable, but this difficulty is not relevant for the present discussion. 

The diagram that I want to calculate is displayed in Fig. 24. I am going to use the mixed 
representation of the propagator, -Dfc(r), that I introduced in previous lectures. I shall call 
the self-energy n(r, p), and I shall focus on the regime where p is a soft momentum while 
the loop integral is dominated by hard momenta /c ~ T. I get ffist 
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FIG. 24: Self-energy diagram in (j)^ scalar field theory 



n(r,p) = j A±.D^^t,{r)Du{-r) (r^ -n^r) (178) 
Dfc(r) = — [(1 + nfc)e-'^'=l^l + n^e'^'^l^l] , (179) 

so that, for r > 

/d^k 1 1 
(27r)^ 2a;p+fc 2bJk 

x[(l + nfc)e-"'=^ + nfce"'="]. (180) 

Now, I take the Fourier transform by integrating over the imaginary time from to /3 

d'^k 1 1 r 1 + np+fc + nfe 1 + np+fc + Uk 







2 /" (i^fc 1 1 J- 1 + rip+fc + nfc ^ 1 + Up+k + ^fc 



^ ^fc — ^p+fc _^ ^p+fc — fik . Qg]^'] 



At this point, I can perform the analytic continuation, iUn — w. In doing so I may run into 
trouble because the denominators may vanish. When the denominators vanish, a priori the 
integral blows up, but this singularity is associated with well understood physics, that I am 
going to discuss. Look at the first term in the integrand 

l+Up+k + rik ^ 1 + rip+k + nk ^^^^^ 



Up+k + (^k- i^n ^p+k + Wfc - (w ± irf) ' 
When the denominator vanishes, it is a signal that there is a process which is allowed. Here 
it is the process by which an excitation carrying momentum p (and energy u) decays into 
a set of two excitations carrying momentum k and momentum p + k. This translates into 
an imaginary part, which is proportional to the rate of such decay. To get this imaginary 
part, we add a little imaginary part to u. In other words, in the continuation, we start from 
Matsubara frequency and we continue up to the real axis but stop a little bit below or above 
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FIG. 25: The various processes contributing to Eq. (181). 

(depending on which propagator we want to consider, retarded, advanced, etc). We just do 
not touch the real axis (see Fig. 23). 

Now let me turn to the numerators, and focus again on the first process in Fig. 25. If 
you view the imaginary part as part of a rate calculation (the denominator providing the 
delta- function that expresses energy conservation), the numerator accounts for the statistical 
factors that accompany the direct process, by which the mode with momentum p decays 
into two other modes, as well as the reverse process by which two excitations recombine 
to form the initial excitation (the second term in the first line of Fig. 25). The direct 
process is enhanced by the factor (1 + nfc)(l + 7^p+fc), while the reverse process will be 
simply proportional to rikUp+k, the probability that the modes k and k + p are occupied. 
Subtracting the two yields 

(1 + nk){l + rip+fe) - nuUp+k = 1 + rifc + ?^p+fc, (183) 

which is indeed the numerator in the first term of Eq. (181). You may verify that all the 
other terms can be obtained from this simple reasoning. 
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Let us consider for instance the third term of Eq. (181). This represents a process by 
which the mode with momentum p can absorb a mode with momentum k and go into a 
state with momentum p + k (see Fig. 25). The reverse process is a mode with p+k going to 
the modes with k and p. The statistical factors here are Uk for the incoming hne {up does 
not count because it corresponds to the particle which I am looking at), and the induced 
emission factor 1 + rip+fc on the outgoing line, giving a factor ?ifc(l + Up^k) for this process. 
There is a factor ?T.p+fc(l + Uk) for the inverse process. You see again that the products 
^k^p+k cancel out, leaving a term linear in n: 

nk{l + rip+k) - rip+kil + rik) = Uk - Up+k- (184) 

Note that this process is a genuine finite temperature effect, and it disappears at T = 0. 
This is in contrast to the the one considered previously, which exists also in the vacuum 
(the 1 in the numerators of Eq. (181)), and which represents a decay process. The processes 
that we consider now are scattering processes involving particles of the heat bath. As you 
will see in a moment, these are the dominant contributions at high temperature. 

Now comes the relations with the rest of the lectures. What we have done so far is an 
exact one loop calculation. Let us pursue a little bit the analysis of the diagram here. The 
loop integral, for the same reason as I discussed already several times, is dominated by the 
largest momenta, that is by momenta of the order of the temperature. Observe the first two 
terms in the self energy, which contain the vacuum contribution, and which represent decay 
processes: 

1 + np+k + nk _^ 1 + np+k + rik ^^g^^ 



^p+k + ^k — i^n ^p+k + Wfc + iUJn 

The denominators are of the order oi k T {p gT), which is big. In contrast, in the last 
two terms, 

^k — npj^k '^p+fc ~ (186) 



^p+k — ^k + 'i^n ^k — ^p+k + 'i^: 

you have difference between two large energies Wp+fc and ujk-, and 

ujp+k-(^k^p- =P-Vk, (187) 

where Vk is the velocity, whose modulus is the speed of light (assuming massless particles). 
Thus the denominator, if u is of order gT, is a soft energy denominator. And indeed, at 
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high temperature, the dominant contribution is obtained form these last two terms 



J [271)^ {2uJk)^ ^ uj + V ■ p —v-p + uj 

r. f d^k 1 r 1 1 ^ dnu 

- 9 / (0 \D(o 7 T^' 188 

J (27r)^ (2a;fc)'' lu — v ■ p cu + v-pJocUk 

where I have used the relation 

'^p+fc-^^fc ~P- ^'^fc = p-'»^fe^^- (189) 

The energy denominators reflect the well-known phenomenon of Landau damping, that takes 
place when the phase velocity of the soft mode, u/p, equals the velocity of the hard particle, 
Vk, in the direction of the propagation of the soft mode. 

What we will do now in the rest of the lecture is to recover similar expressions for QED, 
starting from kinetic theory. My way of showing you the connection with the kinetic theory 
is not a formal way. There is a formal route to deduce things but I shall only be able here 
to give you hints of how things work and are tied together. 

D. Calculation oi IIqet){uj,v) using kinetic theory 

What I shall do now is to do a similar calculation, but for electrodynamics, and using 
kinetic theory. I shall show you how to get the self-energy, Ilq^D{u,p) by solving a simple 
kinetic equation. Doing so, we shall in fact get immediately the hard thermal loop approx- 
imation for IIqed- What is kinetic theory? It is a theory which describes the evolution of 
distribution functions, fg{p,X), which are the phase space densities of particles at space 
time point X = (t, x), carrying momentum p, and electric charge qe (with q = ±1). 

1. Linearized Vlasov equation 

The Vlasov equation is the simplest of kinetic equations. It describes the evolution of 
particles under the action of a force, and reads 

df 

-^ + V-d^f, + Force ■ dpfg = 0. (190) 

The force, in the case of QED, is 

Force = q{E + v AB), (191) 
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where E is an electric field and S is a magnetic field and v = 

I assume that the system is initially in thermal equilibrium, with a distribution function 
which is independent of X: dxP = 0, and function only of the energy of the particle. 

At time t = 0, the system is weakly perturbed away from its equilibrium state, and the 

distribution function becomes 

f-^f + ^f. (192) 

I'm going to assume that the perturbation is small so that I can treat 5f as a small quantity 
and linearize the Vlasov equation in order to determine 5f: 

V ■ dx6fg{p, X) = -qv ■ E^, (193) 



where a covariant notation 



is used. To obtain Eq. (193), I have just replaced in Eq. (190) / by /q + 6f and took 
advantage of the fact that /o is an equilibrium distribution. In particular, I have used the 
fact that /o is isotropic to eliminate the contribution from the magnetic field. 



2. Induced current 
Let me now introduce a new function W{X, v) 

df° 

6f^{v,X) = -qW{X,v)j-. (195) 

If you compare (192) and this expression (195), you can write the following 

f,{v,X) = f^(e,)-qW{X,v)J-. (196) 

You see this is just the beginning of Taylor expansion of the quantity /°(ep — qW{X, v)): 

f,iv,X) = f^iep-qW{X,v)). (197) 

In other words, the distribution function in the presence of the linear perturbation is just 
the thermal distribution function for an energy which is shifted by an amount qW which 
depends on both the velocity of the particle and the coordinates and time. 



76 



The equation in terms of W reads simply 

vdxW{X,v) = v- E. (198) 

This is a first order partial differential equation, which can be solved by the method of 
characteristics. One gets 

W{X,v)= I dt'v- E{x-v{t-t'),t'). (199) 

J —oo 

You can verify by a direct calculation that this satisfies the equation. The interpretation 
is simple. The characteristic line is a straight line which is represented by x = v{t — t'). 
Along the characteristic line the electric field does the work v ■ E in time dt. The work of 
the electric field adds up to make the quantity W in Eq. (199). 

Once we know W{X,v), we can calculate the induced current jj^d(^), which is the 
current generated by the perturbation (in equilibrium the current vanishes). This is given 
by 

jUX) = e I [U{P,X) - /_(p,X)] , (200) 

where the positive charges going in one direction contribute as the negative charges going 
in the opposite direction. By using Eq. (199), we get 

where t = t — t' (not the imaginary time !). 

Now, there is a simple relation between the (retarded) polarization tensor and the induced 
current: 



J'md 



I d'yUZA^-y)My)- (202) 



Since 



dA 

E = -VAo - (203) 

^Ret '^^^ easily obtained from Eq. (201). It is in fact convenient to perform first a Fourier 
transform (F.T.), using 



oo 



F.T. / e-^V(^ - vt) = -4^, (204) 
V - Q + irj 
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(a) (b) 

FIG. 26: Diagrammatic representation of the (one-loop) polarization tensor (left) and induced 
current (right). 

where f{Q) is the Fourier transform of f{X). Then one gets 

nr M ^ ral + . J '^^j^^f^} . (205) 

where 

^1 = ! dpp^^. (206) 

The structure of this equation is very similar to that of Eq. (188). The energy de- 
nominators in particular are identical and reflect the Landau damping processes, while the 
integration over the modulus of the hard momentum (which factorizes) involves the deriva- 
tive of the equilibrium distribution function. The Debye mass (206) is the analog of the 
thermal mass of the scalar field. 

Let me summarize what I have done. I have obtained an approximate expression for the 
QED polarization tensor which could have been obtained also by calculating the one-loop 
diagram in Fig. 26(a), and doing the appropriate kinematical approximation valid when the 
external momentum is soft. Instead of calculating approximately a diagram, I have used 
a kinetic equation to calculate the induce current, from which the polarization tensor was 
obtained. In the kinetic theory, the hard particles are described by a distribution function 
f{p,X) whose slow variations in space time describe the soft, collective, excitations. You 
see how kinetic theory manages to separate the hard and soft degrees of freedom: the 
hard degrees of freedom are those which govern the momentum dependence of f{p,X). In 
f{p,X), p is a hard momentum. The slow degrees of freedom are described by the slow 
variations of the distribution function in space and time (i.e., the dependence of f{p,X) on 
X). In a way, the kinetic theory based on the Vlasov equation can be viewed as an effective 
theory for describing real time phenomena, somewhat analogous to the (Euclidean) effective 
theory used earlier. 
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As I said before, it is possible to establish the correspondence between the kinetic theory 
and the diagrammatic hard thermal loop calculation. Such an approach has been general- 
ized to QCD, where it provides a microscopic effective theory for the quark-gluon plasma 
in the regime where the separation between the hard and soft degrees of freedom can be 
meaningfully realized (see [13] for more details). 



79 



V. LECTURE V 
Introduction 

Today we start exploring a new system, the dilute Bose gas. The specific phenomenon 
that I want to discuss concerns the effect of weak repulsive interactions on the critical tem- 
perature. As you know, in a Bose gas at high density or low temperature, there is a phase 
transition, called the Bose-Einstein condensation. This phase transition occurs in the ab- 
sence of any interaction among the atoms. The question I shall address is what happens to 
this phase transition when the atoms repel each other very weakly. This problem is inter- 
esting in many respects. In the particular context of these lectures it is interesting because, 
as you will see, the techniques employed to calculate the shift of the critical temperature are 
very similar to the techniques that I explained to you in the previous lectures when we dealt 
with general aspects of quantum fields at finite temperature. In particular the technique of 
dimensional reduction, and of effective field theory will play an important role. 

Let me be now more specific. The interaction between the atoms will be characterized 
by a scattering length. I am going to assume that atoms interact only in the s-wave partial 
wave, and I will denote the corresponding scattering length by a. This has the dimension of 
a length. There is another parameter with the dimension of a length, namely n~^/^, where n 
is the density of particles. I am assuming that the gas is dilute, which means that somehow 
n is small. To be more precise, let me remark that together with a and ra, I can form a 
dimensionless parameter, n^^^a. Then the diluteness condition reads n^^^a <^ 1. 

Two issues arise in addressing the question that I was mentioning a minute ago. The 
first issue is whether the phase transition which is observed in the absence of any interaction 
survives in the interacting system. I am simply going to assume that this is the case. The 
second issue, assuming that the phase transition indeed takes place, is connected to the 
evaluation in the change of the critical temperature, and this is what I want to calculate, 
that is 

AT, _ (T, - r°) 
Tc TO ' 

where T, is the critical temperature in the presence of the interaction, while T^ is the 
corresponding critical temperature (i.e., at the same density) in the absence of interaction. 
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What I will show you is that 

AT 

^ = c(anV3), (207) 

c 

where c is a dimensionless positive number. 

This is a non trivial result. The shift in is a quantity which is small if the scattering 
length is small, which is the case when the interactions are weak. Now, given that AT^ — ?■ 
as a — )■ 0, you could think naively that the change of the critical temperature can be 
calculated by perturbation theory. However perturbation theory is useless - well, I mean, 
not completely- but strict perturbation theory is useless because if you start calculating the 
Feynman diagrams order by order in an expansion in powers of a, you will meet infrared 
divergences. That is a situation we have already met. It is an indication that we are doing 
something wrong. What I shall do in the next couple of lectures is to show you the origin 
of the difficulty, and present the techniques that can be used to overcome it, namely the 
techniques based on effective field theory. The exact renormalization group sheds a more 
complete light on this problem, but I shall not have time to discuss it. 

A. Bose-Einstein Condensation 

1. Non-interacting uniform systems 

Let me now remind you of a few basic facts about Bose-Einstein condensation. This is 
textbook material, so I shall skip many details. I consider a collection of non-relativistic 
atoms. These atoms are (spinless) bosons, and at finite temperature the average occupation 
of the single particle level of momentum p is given by 

= e(-p-i/^ - r '- = h 

where fi is the chemical potential. 

Let us start by considering a dilute gas at high temperature. Then the chemical potential 
is negative and large, so that the factor e'-'^p"^-'/"^ is large, and I can approximate rip a; 
g-(ep-M)/r^ In this regime, the gas is essentially classical, the effects of quantum statistics 
(the —1 in the denominator) can be ignored. The density is easily obtained 

n= I l^^p^e"^''^"'' (209) 
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where 

A - (210) 

is the thermal wavelength, an important length scale in the problem. (As in most of these 
lectures, I am using the natural units with h = 1.) It is convenient to rewrite Eq. (209) as 
a formula for fi/T (valid when nX^ ^1): 

M = _in(nA^). (211) 

Let us now examine what happens when one decreases the temperature, keeping the density 
fixed. The formula above gives you the trend (when the gas is very dilute so that nX^ ^ 1): 
the chemical potential decreases in absolute value. As the temperature continues to decrease, 
it eventually reaches the value /i = 0. When /i = 0, the number of particles is given by 
(quantum statistic can no longer be ignored then) 

"-/|^?;7^-^- C(3/2)^ 2,612, (212) 

At that point, something happens. The chemical potential can no longer increase: if it 
would, the statistical factor would become singular (with in particular negative values for 
ep < /i). What happens is that particles start to accumulate in the state with vanishing 
momentum p = 0. Thus, the total number of particles is split into two contributions: one 
contribution, no{T), from particles in the state p = 0, another contribution from the particles 
populating all other momentum states. That is 

n = no(T) + (213) 

Equation (212) can be viewed as the condensation condition relating the critical temperature 
T° to the density n^. This leads to the phase diagram drawn in Fig. 27, where the phase 
boundary is the curve that relates the critical density to the critical temperature T°: 
because A ~ T'^^"^, we get from Eq. (212) oc (T^)^/^. Above this phase boundary, we 
have the condensed phase and below we have the normal phase. It is not difficult to verify 
that the density of particles in the condensate is given by 

r / T \ 3/2i 

no(T) = nl- — , (T<r°), (214) 



TO / 

c 

where n is the total density. The condensate density vanishes, as it should, at the critical 
point. At zero temperature all particles are in the condensate. 
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FIG. 27: Phase Diagram of Ideal Bose Gas 

At this point I want to empasize that the phase transition that we have just described, and 
which occurs in the absence of interaction, exhibits several unphysical features, which will be 
cured by interactions, however small these may be. Consider for instance the compressibility 
(at fixed temperature). This is given by 

1 dV 



^ VdP 



' *^-np(l + np). (215) 



This integral is dominated by the low momentum region, where I can replace Up in the 
vicinity of = by 

2mT° 

p->0. Up ^ (/i = 0), (216) 

an approximation which I have already used in previous lectures. When p — )■ 0, the integral 
is infrared divergent. So the compressibility diverges at the transition, which reflects the 
existence of anomalously large fluctuations of the density. As we shall see, this pathological 
behavior will disappear in the interacting system. 



2. A first look at the effect of interactions 

As I already mentioned, I am working in a system where the interaction is dominated by 
s-wave scattering. I am going to assume that this interaction can be described by an effective 
two body potential, function only of the distance between the two atoms, and of the form 
V{ri — 7^2) = g5^^\ri — 7^2). The relation between the coupling strength g and the scattering 
length a can be obtained by solving the scattering problem and is g = ATia/m. (In fact, 
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FIG. 28: Self-energy diagram of leading order in a 

there are ultraviolet divergences in the scattering calculation, so an ultraviolet cutoff needs 
to be introduced. The coupling strength should be considered as a function of this cut-off so 
that the relation to the scattering length remains valid for any choice of the cut-off. I shall 
not discuss this in detail here since this play no role in our main discussion. I just want to 
alert you about this subtlety.) The hamiltonian density is composed of the kinetic energy 
term and the interaction term 

H{r) = -1-Vtjj^{r) ■ VV^(r) + ^i,^{r)^^{r)^{r)ij{r), 9 = ^, (217) 

where ip{r) and ■?/;''" (r) are quantum fields which describe the atoms and obey the usual 
commutation relation [ip(r),ilj'^(r')] = S^^\r — r'). 

As I indicated at the beginning of the lecture, we expect the shift in the critical tempera- 
ture to be proportional to the strength of the interaction, i.e., proportional to the scattering 
length a. It is therefore natural to try and estimate it using perturbation theory at leading 
order in a. Let us then calculate the correction to the single particle energy /2m due to 
the interaction. This is obtained from the simple Feynman diagram shown in Fig. 28. Let 
me just give you the result: 2gn (it is proportional to g, and the loop integral is proportional 
to the density, so the only hard work is the determination of the factor 2). In the presence 
of the interaction, the single particle energy becomes ep = + '^gn, where = /2m. The 
main feature of this correction is that it is a constant shift, by the quantity 2gn. What is the 
effect of this shift on the transition temperature? Does it move the transition temperature 
up or down when the interaction is repulsive, i.e., (7 > 0? That is, is the shift positive, neg- 
ative, or zero? The right answer is zero! Let me explain to you why this is zero. Remember 
that the criterion for condensation is n(yU = 0,T) = 2.612, where 
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This is for a = 0. When a 7^ 0, the density is given by 



n = 



I 



£p 1 

(27r)3e(^p-(/^-29"))/'r - 1 



n{fi — Afi,T), Afi = 2gn, 



(219) 



where fi — A/i = fi — 2gn can be viewed as a modified chemical potential since 2gn is a 
constant. The condensation condition that we had for a = 0, n{0,T) = 2.612, is unaffected; 
it just occurs at a different chemical potential (/i = 2gn), but the relation between 
the critical temperature and the critical density is independent of g. So the first-order 
perturbation theory or, if you wish, the mean-field calculation, does not produce any shift 
in the critical temperature. This result, simple to establish, shows immediately that the 
linear relation between the shift in Tc and a is not going to be obtained by a trivial procedure. 

Q: If you have temperature dependence in the mean- field, somehow... 
A: Somehow? What do you mean by "somehow"? 

Q: If you calculate only this diagram, of course you do not get a temperature dependence. 

A: You can think of a more complicated calculation. But at this level, the message is fairly 
robust, although it has in fact been overlooked by quite a number of people. Perhaps 
I should say that I am approaching the transition temperature from above where the 
gas is classical, and there is no... Well let me not say that now. I was about to say 
"symmetry breaking", but I have not mentioned that concept. I shall do that shortly. 

Although in leading order the transition temperature is not affected by interaction, the 
interaction has nevertheless a profound effect on the properties of the system. Let me show 
that in the case of the compressibility. I leave it as an exercise to you to recalculate 
Let me just give you the result. 



(You see why you get the denominator: when you take the derivative with respect to fi of 
the density 



1 dn 5^L^p(l +^p) T 



(220) 



n^dfi I + ^ J^np{l + Up) 2g 
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you have an explicit derivative, but you have also to take the derivative of the factor n inside 
the integral, and that brings another factor dn/dfi (proportional to 2g/T).) The momentum 
integral in Eq. (220) are infrared divergent, but now the divergences actually cancel out since 
the same divergent integral appears in the numerator and in the denominator. Of course, if 
g is very small, the compressibility is very big but even for infinitesimal g the compressibility 
becomes finite. So there is a deep modification of the properties of the system. What is being 
done here, in terms of Feynman diagrams, is actually a resummation of a chain of bubble 
diagrams which contribute to the screening of the long wavelength density fluctuations. 



There is another aspect which is qualitatively new in the presence of the interaction: we 
can discuss the phase transition in terms of symmetry breaking. I am going to consider very 
low, in fact zero, temperature and show you that Bose-Einstein condensation is what one 
may call a quantum phase transition seen here as a change in the properties of the Bose 
gas as one tunes the chemical potential. To carry out the discussion in simple terms, let me 
introduce a quantum state, which I call |iVo), 



This is a coherent state, with aj the zero momentum component of the field creation operator, 
[ao,aQ] = 1. The expectation value {NqIH — /iaJaojA'o) contains no contribution from the 
kinetic energy because carries no momentum. The result is 



If g = and yU < 0, the ground state in Fock space is A'^q = 0: this is the vacuum state 
with no particle. On the other hand if = the system is completely degenerate: we 
have an arbitrary number of particles in the ground state because it cost no energy to add 
one. This degeneracy is the source of the large density fluctuations in the non interacting 
system at the transition. If > 0, the system is completely unstable: you can decrease 
its energy by an arbitrary amount by adding more and more particles. The presence of an 
interaction controls this phenomenon and cures the instability. When /i > 0, the quadratic 
term coming from the interaction generates a minimum at the value A'^o = fi/g (see Fig. 29). 



3. Symmetry breaking 




(222) 



{No\H - fialao\No) = -/iA^o + |iVo- 



(223) 
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FIG. 29: Energy of the coherent state \No) as a function of A'^o for fi < and g = (straight Hne) 
and for /i > 0, in which case a minimum occurs for A'o = n/g- 

condensation vacuum 

1 ^/i 

No = M/g 

FIG. 30: Phase diagram of the weakly Interacting Bose Gas at zero- Temperature 

As a function of fi, the phase diagram looks Uke Fig. 30. This is a pattern which allows us to 
view this Bose-Einstein condensation using techniques which are familiar in quantum field 
theory, namely borrow all what we know about symmetry breaking and associated Goldstone 
modes and related phenomena. The symmetry which is broken here is a U{1) symmetry: 
Remember that the hamiltonian of the system does not depend on the (global) phase of the 
bosonic field. However, the field acquires an expectation value (that is explicitly taken into 
account here in our description with a coherent state), the U{1) symmetry is broken. 

Q: Excuse me, did you include —fiN to hamiltonian? 

A: Yes. Whenever I talk about the hamiltonian H, it is the true hamiltonian with kinetic 
energy, potential energy and —fiN. When I talk about degenerate states, I am referring 
to eigenstates oi H — fiN. 

This is a very rapid digression, but it indicates to you that indeed even an infinitesimal 
interaction changes the properties of the system in a qualitative way (you can extend this 
discussion to finite temperature). In the presence of interactions the ground state is indeed 
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different and we can use languages which we are famihar with in field theory to understand 
what is going on. This actually leads us to deep consequences; when there is a symmetry 
breaking there are Goldstone modes associated with a rotation of the phase of the order 
parameter and those Goldstone modes have a strong impact on the physics of Bose-Einstein 
condensates. But I shall not discuss this too much here. 

There is another digression that I want to make to illustrate the role of mean-field effects 
on Tf.: these are indeed very different in uniform and finite sytems. 

4- Atoms in a trap 

I would like to discuss briefly what happens for the atoms in a trap. This for two reasons. 
The first is that the physics of cold atoms in traps is what has triggered the renewal of interest 
in Bose-Einstein condensation and much of the works that I am discussing. As you know, 
one is able now to cool atoms in a trap to a sufficiently low temperature to observe Bose- 
Einstein condensation, in a system where the interaction strength is very small. Till then, 
the prototype of systems in which one could observe a phase transition akin to Bose-Einstein 
condensation was liquid helium. But in liquid helium the interaction between the atoms is 
very strong and the fraction of the particles that are sitting in the condensate is never bigger 
than about 10 %, even at very low temperatures. Nowdays, one is able to prepare small 
condensates of few tens of thousands of atoms in a trap, tune their mutual interaction to 
be as small as desired, and observe genuine Bose-Einstein condensation. However, and this 
is the second motivation for this digression, the presence of the trap inhibits some of the 
effects of the interactions that occur in uniform system. So the study of atoms in a trap 
gives, so to speak by contrast, an interesting perspective on some aspects of the effects of 
the interactions, related to long wavelength phenomena that are characteristics of uniform 
systems. These are these long wavelength phenomena that, in my view, make Bose-Einstein 
condensation so interesting. 

I am going to consider the following situation: the trapping potential is described by a 
harmonic oscillator, with a typical level spacing hw. I assume that the temperature T is 
high, such that k^T ^ hu. Under this condition, I can use a semiclassical approximation to 
describe the atoms. What does it amount to? The semiclassical approximation essentially 
states that one can consider the gas of the atoms locally as a piece of uniform matter whose 
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density is equal to the local density in the trap. In other words, one assumes that the energy 
of an atom can be written as 

ignoring the fact that p^ and do not commute. The density is, as usual, obtained by 
integrating the distribution function over the momentum of the particles 

n{r,n,T) = / f ^ \ — -. (225) 

J (27r)3 e(^(P''^)~^)/^ — 1 

This is now a function of r. The total number of particles is 

N = J £rn{r,fj.), (226) 

and is kept fixed. The number density at the center of the trap, where it is is the biggest, 
is given by 

"^°'^'^) = / (^ e(-(^.o)-LT_i - (227) 

When 

n(0,/i = 0,T) = 2.612, (228) 

the atoms at the center of the trap will undergo Bose-Einstein condensation. That is, 
condensation occurs whenever the density at the center of the trap satisfies the same relation 
with temperature as in a uniform system. The density at the center of the harmonic trap 
can be calculated easily. You find then, that in the absence of interactions. 



(229) 



huj VC(3), 

So you see that if is big enough (in typical experiments, ~ 10^ ^ 10^), the condition 
for the validity of the semiclassical approximation, ksT ^ hw is well fulfilled. 

Now we can discuss easily the effect of the interaction on the transition temperature of 
the trapped gas. I can ask again my question on the shift in T^. Will it be positive, negative 
or zero? Let me tell you that this is not going to be the same result as before, namely, 
this is not going to be zero. That leaves two possibilities: positive or negative. You should 
be able to answer this question because that requires no calculation at all. What will the 
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interaction do if you put particles in the trap? Let us imagine putting particles in a trap, 
and keep adding particles (at fixed temperature) until the density at the center of the trap 
is high enough for condensation to occur. When the particles in the middle of the trap just 
begin to condense, let us switch on the interaction. What happens? Because the particles 
repel each other, the gas will expand and the density in the middle of the trap will decrease, 
destroying the condensation. How can one recover the condensation? By decreasing the 
temperature: this will indeed decrease the kinetic energy of the atoms making them more 
sensitive to the effects of the trapping potential (that pushes them towards the center of the 
trap). This discussion shows that the effect of the repulsion between the atoms is to shift 
the transition temperature downwards. This effect has been observed in experiments. It's 
magnitude has been estimated 

1.32— iV^/^ (230) 



Tc ah, 



o 



where Oho = ^/h/mcu. 

This negative shift of Tc is purely a mean-field effect, and its physical origin is transparent, 
as we have seen. 



B. Towards the calculation of ATc in an uniform system 

In a uniform system, the effect of the interaction leads to a positive shift in the critical 
temperature, which is therefore opposite to the effect that we just discussed for atoms in 
a trap. The physics responsible for a change in Tc in uniform systems has in fact nothing 
to do with what happens in a trap where mean field effects dominate: as we have seen, in 
uniform systems mean field corrections do not produce any shift in Tc. 

Let me first of all indicate a useful relation, which is a purely geometrical relation, valid 
in leading order in the interaction strength. It is a relation between the shift in the critical 
density and that in the critical temperature. Figure 31 represents the phase diagram in the 
density-temperature plane in the cases where a = and a > 0. I am assuming here that 
in the presence of the interaction, there is a still a phase transition similar to Bose-Einstein 
condensation, and that if a is small the critical line is only slightly displaced from where it 
is when a = 0. The relation that I am referring to is 

= -, 231 

Tc 3 ric ' ^ ^ 
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FIG. 31: The critical lines nc{Tc) for the non interacting (n[?) and the interacting (a > ) systems. 



which just follows from the fact that the two curves are very close to each other, and the 
fact that I am looking at the leading order in a. The factor 2/3 has its origin in the 
relation n° ~ (T°)'^/^, while the minus sign is obvious from the figure. Why is that relation 
important? It is important because Aric is much easier to calculate than ATc. This is so 
because it is easier to calculate at a fixed temperature than to calculate at a fixed density 
(in the latter case, you need to adjust the chemical potential to keep the density fixed as 
you change a). 

Let us go back to the interaction. The hamiltonian that we consider is of the form 

H = j d'r l^t(^) (^_^ ^^r) + |v;t(r)^t(^)^(^)^(^)| , 9 = ^, (232) 

which is the hamiltonian (217), and we have recalled the relation between the coupling 
constant g and the scattering length a. When using this hamiltonian we assume that we can 
ignore a lot of details of the atomic physics. The potential between the atoms is replaced 
by a contact potential, which is of course meaningful only if the atoms are on the average 
far from each other. Only then can we ignore the details of the atom-atom interaction, its 
dependence on specific electronic levels, etc. To be more precise, this hamiltonian will be 
used to describe modes of the Bosonic fields whose wavelengths 1/k are large compared to 
the range ro of the potential. 



Zero < 1. 



(233) 



91 



We have also mentioned another condition, that involves the scattering length a, 

na^l-' < 1. (234) 

This is the statement that the distance between the atoms is large compared to a. Recall 
finally that in the vicinity of the transition, riA^ ~ 1, where A is the thermal wavelength. 
The condition on the scattering length translates then into the relation a ^ A. 

In order to calculate Aric, it is useful to express the density in terms of the propagator 
(in the imaginary time formalism): 

n = lim {T^p{T,r)ij\0,r)) = lim G{r,r). (235) 

The Fourier transform of the propagator obeys the Dyson equation: 

G^\iUn,p) = GQ\iujn,p) + T.{iuJn,p), (236) 

where Go is the free propagator, S the self-energy, and Un = 2miT is a Matsubara frequency. 
Therefore, the expression for the density of particles can be written as 



-ibJnT 



n = limTy; / , (237) 

n 



r->o- {2n)^ ep - icun + ^{iUJn,P)' 

where ep = — /i and Gq ^(ia;„,p) = 6° — /i — iun- This is an exact relationship: It allows 
me to calculate n, provided I know how to calculate the self-energy. 

Now I need to say a few words about the condensation condition. I shall give you the 
result and try to motivate it. I am going to assume that condensation takes place when 

G^\uj = 0,p = 0) = 0. (238) 

Note first that, in the non interacting case, the condition Gq ^(0, 0) = yields /i = 0, which 
is indeed the condensation condition that we have already met. Turning to the interacting 
system, one may recognize that G~^{iUn = 0,p = 0) is just the second derivative of the 
"effective potential" (the free energy expressed in terms of the expectation values of the 
field ip and tp"^). Assuming that the condensation is a second-order phase transition, the 
second derivative of this effective potential vanishes at the transition. This is just the 
statement that G~^(0,0) = 0. 

By using the explicit expression G~^{iuJn, p) = e^ — fi — iUn + S(za;„, p), one sees that the 
condition G~^ = reduces to 

/i = S(0,0). (239) 
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FIG. 32: Second order diagram contributing tho the self-energy. 



Now I can calculate the shift Arte in the critical density. I get 

An = lim rVe"*"'"^ / _fp_ ( 1 1 1 .240) 

' r^o- ^ y (27r)3 \eO + S(zw,p) -S(0,0) e° - J ' ^ ^ 

The first term is the critical density of the interacting system, for which fi = E(0,0), the 
second term is the critical density of the non interacting system at the same temperature, 
for which fi = 0. Note that the only place where the interaction enters is the self-energy. So 
we have to evaluate the self-energy. You see here that when the self-energy is frequency and 
momentum independent, there is no correction, because if the self-energy does not depend 
on or p it is the same as for u = p = 0, and then the two terms cancel and Auc = 0. 

Q: This is precisely what happens in the mean- field approximation? 

A: Yes. As I said, the mean field approximation leads to a self-energy which is independent 
of frequency and momentum. That produces no shift of T^. To get a non-trivial effect 
you need frequency and momentum dependences in the self-energy. 

Let me start the calculation of the correction. We have already done the first-order 
calculation and shown this to be zero. Therefore, we have to go to the second-order. I shall 
come back to it next time. The second-order diagram is shown in Fig. 32. The corresponding 
contribution to the self-energy reads 

731, /• I I I 

(27r)3 ek^q - i{Un' - efc - iUn' ep+q -i{uJn + ' 

(241) 
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where the factor comes from the double sum over the Matsubara frequencies. I am going 
to focus on the contributions with n = n" = 0. Remember that the infrared divergences 
of the scalar field theory are coming from this particular sector where all the Matsubara 
frequencies vanish. I am going to look at this particular contribution, and follow the same 
strategy as a few lectures ago when I discussed the divergences in high temperature QCD 
and I exhibited special classes of Feynman diagrams which are infrared divergent. I used a 
power-counting argument and replaced the multiple integrations over the momenta by an 
integration over a big momentum vector in a large space and I ignored the angular integral. 
Let me proceed in the same way quickly today (I shall come back to this result next time) 
because I would like to show you that there is difficulty. Think therefore of the momentum 
integrals over q and k as an integral of a big vector in 6 dimensions. Then the integral in 
Eq. (241) behaves as 



This is a logarithmically divergent integral. 

Q: How about the external momentum p? It does not provide any cutoff? 

A: Well, it does, but in a subtle way that we shall discuss more precisely next time. For 
today, let me just observe that what we need to calculate is S(0,p) — S(0,0), and 
there is no p in S(0, 0). 

What will come after is very much the same thing as what we met when we discussed 
QCD at finite temperature. Namely, when I calculate the third-order, the fourth-order, etc., 
I find increasingly divergent contributions. The pattern of divergences becomes worse and 
worse as the order increases. What I want to argue next time is that this situation is very 
much similar to what we have met already. You will see that this pattern of divergences is 
indeed that of 3-dimensional field theory. Having recognized that, we will know what to do, 
namely, how to construct an effective theory for the sector where the divergences occur. This 
will allow us to establish the linear relation between AT^ and a, and to obtain an explicit 
formula for the coefficient c in the formula (207). I will then explain to you why it is hard 
to calculate explicitly this coefficient. 

Q: In the last integral, I think you used a large momentum approximation? 




(242) 
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A: Not really. But, of course, the logarithmic divergent integral involves a ratio of momenta, 
and it is large when one of the momenta is large (or small) compared to the other. 
But I am really concerned here by the low momentum sector. 

Q: This approximation is on the low momentum? 

A: Yes. I am assuming that all the loop momenta are going to zero at the same rate. Note 
that, for this second order case, I can calculate the integral explicitly and I shall do 
that next time. I hope then that this question will be completely clarified. 
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FIG. 33: The critical lines nc{Tc) for the non interacting (n[!) and the interacting (a > ) systems. 

VI. LECTURE VI 

Let me remind you that we want to establish the following formula for the shift in the 
critical temperature of the Bose-Einstein condensation caused by weak repulsive interactions: 

AT 

^ = c{an'/'), an'/' « 1. (243) 

c 

We want to establish this formula to leading order in the strength of the interaction, mea- 
sured by the scattering length a. The plot in Fig. 33 displays the critical line for the 

3/2 

non-interacting gas. It goes like nc ^ Tc . Also drawn is the critical line for the interacting 
system. We assume that it differs very little from that of the non-interacting system when 
a is small, and the curve is drawn here for a positive scattering length. 

I'm interested in the critical temperature Tc at point C on Fig. 33: Tc is the transitiion 
temperature of the interacting system for some small positive value of a. Starting from the 
non interacting system, a = 0, I can approach this point C either by moving the density 
downward at constant temperature (going from A — )■ C) or moving upward the temperature 
at constant density (going from B ^ C). I argued last time, based solely on the geometrical 
propeties of Fig. 33, that AT^, and Anc are related by the simple equation: 

ATc 2 Aric 



Tc 3 nc 



(244) 



I also argued last time that it is much easier to calculate Att-c, i.e., to work at fixed temper- 
ature. 

The last thing I want to remind you is the condensation condition, expressed as the 
vanishing of the inverse propagator for vanishing momentum and Matsubara frequency. 
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This translate into a condition on the self-energy: 

^"^(a; = 0,p = 0) = ^ S(a; = 0,p = 0) = (245) 
With this, we can write the formula for the shift Aric: 

(2vr)3 

/ 1 1 

\eO + S(2a;„,p) - S(0,0) - - m 

The first term of the integral is the density of the interacting system at criticality. The 
subtracted term is the density of the same system at the same temperature, but without 
interaction (we used here the fact that at criticality the chemical potential in the non- 
interacting system is zero). From this formula you may observe that if the self-energy is a 
constant, independent of frequency and momentum then S(2a;„,p) — S(0,0) vanishes, and 
so does Aric- This is what happens for instance in the mean field (one-loop) calculation that 
I did last time. 

Now perhaps I should add some word of caution here. Sometimes you may be led to 
consider a mean field calculation where the dominant effect of the interaction leads to a 
momentum dependent potential, that is, S doesn't depend on frequency but depends on 
momentum (note that for this to occur, you need to go beyond the approximation that 
consists in replacing the atom-atom interaction by a contact potential). That dependence 
on momentum can often be described by an effective mass. Such an effect would produce a 
shift, just because the effective mass enters for instance the thermal wavelength. But this 
is not the kind of effect that I am interested in here. 




A. Power counting and infrared divergences 

We also started discussing the interaction. I remind you that the typical hamiltonian we 
use is given as 

H = Jdh jv^^r) (^-^^ ^(r) + |v^t(r)^t(r)^(r)^(r)| , 9 = ^- (247) 

We argued that the leading order correction to the single particle energies, of order g, 
produces no shift in Tc, and we started considering the calculation of the second order 
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FIG. 34: Second order diagram contributing to the self-energy. 

diagram (see Fig. 34). Then I argued that one could expect infrared divergences when all 
Matsubara frequencies are equal to zero. So, I'm going to calculate this diagram only in this 
particular case. Let me give you the complete expression, near criticality: 

d^A; f d^g 



Eitoon = 0, p) - S(0, 0) = -2/r2 J J 



(2vr) 



3 



where /i' is given as 



1 /I 1 



^i' = fi-2gn = -^ (249) 
2m 



I'm doing a calculation in second order perturbation theory. But I'm incorporating in 
the chemical potential the shift due to the first order diagram. At the mean field level, 
the condensation takes place in exactly the same way as in the non-interacting system: 
/i' is negative until one reaches condensation where it vanishes. The notation with Kc is 
convenient, as it allows me to express the energy denominators as follows: 

e'^-f,' = ^{e + 4), (250) 

and Kc will play the role of an infrared regulator in the integrals that I'm going to calculate. 
Within the present perturbative setting, it measures the deviation from criticality. 

Let me give you the value of the integral in Eq. (248). This can be calculated analytically 
for finite Hc'- 

.,0.p,-.,0,0)^l^(^)'{5^.^-.^.i,n(l.(£-)^)-| ,51) 

This result allows me to do a couple of comments. The quantity S(p) — S(0) is plotted 
schematically in Fig. 35. For small p, p ^ k^, it goes like p^, while for large p, p ^ Kc, 
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FIG. 35: — S(0) given by Eq. (251) for small p and large p 



it goes like In p. When p is very large compared to Kc, S has a logarithmic behavior. And 
when p is small compared to smooth function which has a regular momentum 

expansion. Of course Hc is a quantity that we would like to let go to 0, because we want to 
approach the condensation. But when goes to we get an infrared divergence. That's 
the infrared divergence that I mentioned at the end of the last lecture. 

Now I will show you that this peculiar infrared behavior is actually not limited to the 
second order diagram. It will show up in higher orders. We are going to rediscover the 
pattern that we have already identified in the case of hot QCD, or in the case of the thermo- 
dynamics of the scalar field that we discussed a few lectures ago. Let me repeat the analysis 
that I did last time for the second order diagram (Fig. 34). Focusing on the region where 
all the momenta are going to zero at the same rate, we regroup the two loop-momenta into 
a single 6-dimensional vector K, and write 



This is a two loop diagram, so there are two momentum integrations, hence the d^K. Then 
are three propagators of the form 1 / [K"^ -|- k^)^. And there are a number of coefficients: there 
is a factor g"^ coming from the two vertices, and a factor coming from the two summations 
over the Matsubara frequencies ( I keep only the zero Matsubara frequencies, but the factor 
that accompanies the double sum remains). The factor (2m) comes from the p^ /2m in 
the three propagators. The last line is obtained by remembering that g ~ a/m, and that 
mT ~ Also, to make it easier to see the divergence of the integral when 0, 

I have removed n^. from the denominators, and put it at the lower end of the integration. 
Then I ignored the angular integration, which would produce just a numerical factor. You 




(252) 
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FIG. 36: 2+1 loop diagram 

see now that in the hmit where Kc goes to this integral is logarithmically divergent. Also, 
it is proportional to (a/A^)^. 

Let us now generalize this calculation. I go from this 2-loop diagram to an /-loop diagram. 
Let's see what happens by adding a simple loop as shown in Fig. 36. When I add a loop I add 
one vertex and two propagators. One can easily verify that this is generic. The calculation 
of the /-loop diagram is then going to take the following form: 



\l-2 



In this calculation I have added / — 2 loops to the 2-loop digram of Fig. 34, in order to get 
an /-loop diagram. By doing so, I have added / — 2 loops, / — 2 vertices, hence the factor 
gi-2rpi-2^ Then I have an integration (d^Ji')'"^ and 2(/ — 2) propagators {K"^ + 
and a factor (2m)^'^'~^\ You are still with me? I can rewrite these factors as 

In the integral, I remove Kc in the denominator, and put it as a lower bound. After all this, 
I can rewrite the expression (253) as 



so that the /-loop contribution to S reads 

As anticipated, we are indeed confronted to a situation which is very reminiscent to one 
that we have already encountered. Namely, as all terms are infrared divergent. I 

will show you later that there are good reasons to choose of order a/A^. But if you do 
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FIG. 37: Some diagrams contributing to the parameters of the effective theory 



that, a/A^ ■ l/nc is of order 1. That means that all the diagrams in perturbation theory will 
be of the same order of magnitude. Perturbation theory breaks down. It cannot be used 
to calculate the shift of the critical temperature. But we know what to do. This is very 
analogous to what we have encountered earlier. And I have shown you in the case of the 
scalar field that we can develop other tools to handle this problem. One of them is effective 
field theory, to which I now turn. 



Infrared divergences occur in Feynman diagrams only when all the Matsubara frequencies 
Un are set equal to 0. Let's then focus on this particular sector and recall what we did for 
the scalar field. The effective theory amounts to retain in the field expansion 



only the component ipoir), which does not depend on the imaginary time. As you know 
from previous lectures the effective theory for ipolr) will be an effective theory for a field in 3 
dimensions. One dimension is lost because I have abandoned the imaginary time dependence. 

If I just replace the field ?/'(r, r) by ipoir), I can perform trivially the integration of the 
imaginary time in the action, and get: 



B. Effective theory 



(257) 




(258) 



where 




/^)^o(r) + |(|V^o(r)^)^ 



(259) 
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This is the leading order, where the effective action is just /3 times the energy of the field 
configuration. This is what I called earlier the classical field approximation. We know how 
to calculate the corrections to this simple approximation. Remember that I described to you 
the systematic procedure to do that. What we have to do is to calculate Feynman diagrams 
where the external lines corresponds to ipoi^) or V'o(^) ^^"^ ^^e internal lines corresponds 
to the modes with non-vanishing Matsubara frequencies. Such diagrams are displayed in 
Fig. 37. The diagram of Fig. 37a is proportional to a, and it is a correction to the chemical 
potential fi. We shall not need to evaluate this correction because the chemical potential 
is eventually adjusted so that the system is critical, or at the phase transition. But there 
are other corrections. Let me first look at the correction to the coupling constant, the 
diagram Fig. 37b. This correction is proportional to a^. Since we are interested in leading 
order calculation, it can be ignored. The correction in Fig. 37c represents a three-particle 
interaction, and corresponds to a term not present in Eq. (259). This is of order a^, and can 
again be ignored. The lesson of this brief analysis, which we have done more extensively in 
the case of the scalar field theory, is that if one is interested in the leading order effect of the 
coupling constant, one can just use, as an effective theory, the classical field approximation, 
whose energy density is given by Eq. (259). 

However, we are a priori somewhat stuck here, because with this effective theory per- 
turbation theory cannot be used even though the coupling g can as small as one wants. 
Let me briefly review that issue, and recall why a priori we should expect a problem with 
perturbation theory. When we do perturbation theory, we are assuming in some way that, 
in Eq. (259), the kinetic energy 2^(V^/'o)^ is big compared to the potential energy g (l^/^oP)^- 
Naively we expect that if g is very small, the potential energy is correspondingly small. But 
this is not always the case. Perturbation theory breaks down precisely at that particular 
scale where all the terms in the effective action are of the same order of magnitude, that is, 
in particular, when (in average) 



Let me repeat here an analysis that we have done in a more general context. Remember 
that we can separate layers of fluctuations corresponding to different momenta or different 
wavelengths. Consider the particular fluctuations of momentum Kc, so that 




(260) 




(261) 



102 



The density in the effective theory here is (|'?/'oP)- This is also given by 

..~(W'>~/(0; ^,,,i._, . (262) 

Close to condensation, I can ignore the chemical potential (or a constant self-energy), and 
replace (ep — /i)/T by /2mT. For long wavelength modes, p^/2mT is small and I can 
expand the exponential factor, and rewrite Eq. (262) as 

n~ / , , ■ 263 

} (27r)3 p2 ^ ) 

(This approximation is in fact in line with the classical field approximation.) This integral 
diverges at large momentum. But it is supposed to represent only the contribution to the 
density of long wavelength modes, with p < k^. With as ultraviolet cut-off, this gives 
n oc niTKc- Let's then compare ^IV^oP with g (I^/'qP)^, or equivalently ^ and g\ipo\'^, with 
(iV^oP) ~ mTnc- I can rewrite Eq. (260) as 

^ ~ gmTKc. (264) 
2m 

Remembering that g a/m and mT ~ one sees that the condition (260) implies that 

«:c ~ ^. (265) 

This is the characteristic scale that we encountered before and which signals the breakdown 
of perturbation theory. We see here how this scale emerges from a simple analysis of the 
effective action. 

Now, if we cannot expand, what can we do? We can still make progress, because we 
know that Eq. (246) is an exact formula for ArZc, and we can deduce from it an equivalently 
"exact" formula for ArZc in the framework of the effective theory. This is very easy to do. 
We just ignore all the terms which have n different from in Eq. (246). This leads to 

^"= = /(5^{4 + E.(p)-E.(0)-y- 

where Sci(p) = S(2a;„ = 0,p). This formula is "exact" in the effective theory, in the sense 
that it does not involve any approximations, beyond those made to arrive at the effective 
theory. 

I'm going now to introduce a new notation. Let me set 

[/(p) =2m(E,Kp)-S,K0)). (267) 
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p 



FIG. 38: General behavior of U{p) as a function of p. 



Then, with a httle algebra, I can rewrite the formula (266) as 




(268) 



This is still an exact formula (within the effective field theory). And if I know the exact 
expression for U{p), then I can conclude something. Now, I have to tell you a few more 
things about the general behavior of U{p). Before, I have done the calculation with an 
infrared regulator: remember that when the momentum was smaller than the regulator, 
U{p) was analytic, ~ p"^. But when one removes the regulator, the behavior is actually 
logarithmic (in leading order), and more generally a power law, + U{p) ~ p^'"^ where 
7], a small number, is the anomalous dimension (a small number of order 0.1). This power 
law behavior is characteristic of the scaling regime of a second order phase transition. The 
main point at this stage is that the integral in Eq. (268), if the function U{p) behaves as I 
indicated, is a finite integral. 

There is more that we can say, namely that U{p) > 0. 1 do not have a complete, 
analytical, proof to offer for this property. One may just observe that it is verified in the 
second order calculation, Eq. (251). The critical fluctuations correct the small momentum 
behavior {p"^ — t- p^"''), but do not alter the fact that U{p) is a growing function of p at 
small p. As for the large momentum behavior, it is correctly given by perturbation theory. 
Detailed calculations confirm that U (p) is an increasing function of p, with the general shape 
displayed in Fig. 38. Now, if U{p) is positive, so is the integral in Eq. (268). Therefore Auc 
is negative, which implies that ATc is positive. This answers the question concerning the 
sign of ATc. 

The last thing we have to do is to show that the shift of Auc is linear in a. This can 
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be done via a simple analysis. Let me first remark that the effective theory which I have 
written here in principle requires an ultra-violet cut-off. You have seen that already from the 
calculation of the density in Eq. (263): if you want to calculate the density with a statistical 
factor of the form you need an ultra-violet cut-off. This is in line with the philosophy 
of the effective theory. Remember what we did in the scalar field theory: we introduced an 
intermediate scale, i.e., a separation scale A, and we integrated the modes above this scale A 
to get an effective theory valid below that scale. Of course it may happen that the effective 
theory contains ultra-violet divergences, and these divergences are in principle compensated 
by the cut-off dependence of the coefficients of the effective theory. Here you don't have to 
worry about such issues because, as I argued, the density enters mostly the correction to the 
chemical potential. Observe also that in the integral of Eq. (266) we are taking the difference 
of two contributions, and this difference is ultra-violet finite. But there are cases where we 
need a cut-off A. What is the scale of this cut-off? To answer that question we need to 
ask ourselves first why we get ultra-violet divergences in the effective theory. The answer 
to that is contained in Eq. (262). The full statistical factor there, l/(e(^p~^')/^ — 1), kills all 
the momenta which are bigger than the temperature. But we are using an approximation 
to the statistical factor that is valid only at small momenta, i.e., for p < y/mT, so that the 
natural cutoff for the effective theory is of order \/mT ~ 1/A. 

Now the effective theory is a 3-dimensional theory, and is "super renormalizable" . What it 
means is that the only ultraviolet divergence is the one that I have discussed. It corresponds 
to a correction to the chemical potential (or to the mass, in the language of field theory). As 
I have argued, such a correction to the chemical potential is innocuous because the chemical 
potential is adjusted to be at criticality. Let us then go through some dimensional analysis. 
As you've seen, there is a natural momentum scale in the problem which is a/A^. So let us 
set p = Xp- (with X a dimension-less variable). In principle, quantities calculated within the 
effective theory depend also on an ultra-violet cut-off A ~ 1/A. Thus I can write 

^(p = a:^,A) = (^) ^(^'X^xJ' ^^^^^ 

where -^j^ ~ a/A if A ~ 1/A. In fact, since the theory is super renormalizable, I can let the 
cut-off A go to infinity, and get a finite result. But the fact that I can let A go to infinity 
does not guarantee that the results will not depend on A, if it is kept finite. Finite cutoff 
corrections will be truly negligible only when a/A is sufficiently small. In that case, I can 
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replace Eq. (269) by 



U{p = X 



a 



2' 



A) ^ (^)V(x), 



(270) 



where (t{x) is a universal function of x (that is, independent of a). 
Then I can rewrite Aric as 

2 a f°° dx xa{x) 



TxX^ Jq X x2 + (t(x)' 



(271) 



where I have extracted the logarithmic integration measure to emphasize that the variations 
of the integrand are best visualized on a logarithmic scale (see Fig. ?? below). The integral 
is just a number, so that An^ is proportional to a. Basically, the dependence in a just 
follows from dimensional considerations, once we have made sure that possible ultraviolet 
cutoff effects play no role (which, as we have argued, requires a/A to be small enough). 

At this point, it is useful to briefly review the steps that led to this result. The first step 
is to recognize that perturbation theory doesn't work. Naively one could think of expanding 
in powers of a, because a is very small quantity. But we have seen that this doesn't work 
because the calculation of Feynman diagrams is plagued with infrared divergences. We 
have exploited the fact that these divergences occur when all the Matsubara frequencies 
are vanishing. Then, by relying on what we did earlier, I explained to you how one can 
deal with all the Feynman diagrams at once, by constructing an effective theory. In this 
particular case, the effective theory is extremely simple because it just amounts to replace in 
the original action the bosonic field, which is a function of three space coordinates and the 
imaginary time, by a field which is independent of time and depends only on three spatial 
coordinates. This leads to a three dimensional effective field theory, which in principal allows 
us to calculate Arig. 

The effective theory makes obvious the reasons for the breakdown of perturbation theory: 
there exists a typical scale of fluctuations at which kinetic energy and potential energy are 
of the same order of magnitude. Therefore, in order to calculate within this effective theory, 
we have to use tools other than perturbation theory, even when the coupling is small. 

But even without doing explicit calculations, we have been able to extract from the 
effective theory the answers to the questions that we are addressing. Namely the fact that 
(a) Aric is negative or ATc is positive and (b) Anc is linear in a. 

Q: You take some kind of classical field approximation somewhere? 
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A: Yes. The effective theory that I have discussed is what I called the classical field ap- 
proximation. I have argued that Feynman diagrams are divergent when all Matsubara 
frequencies are vanishing. Earlier in the lectures I have shown that we can construct 
an effective field theory to handle systematically all these Feynman diagrams at once. 
This effective theory is in general an infinite series in local operators, such as (f)^, (p^, 
0^, 0^(V0)^, etc. An important point is that the coefficients can be calculated in 
perturbation theory. If you are interested in leading order, then the leading terms 
in the expansion are enough, and this leads to the classical 0^ theory in 3 dimension 
(with 0(2) symmetry). So, in the particular context of the present problem, things are 
extremely simple because the effective theory is just the leading order term. Namely 
you just take the initial action, which is 

^ = jy^J r{r, r) ^(r, r) + ■ ■ ■ , (272) 

and replace ?/'(r, r) by 'ipQ{r), to get 

S = P j rfV^*w(-^)^o(r) + ---. (273) 

Q: Is that a classical field approximation? The Heisenberg h goes away by this procedure. 

A: Yes. If you wish. Another way to think about this approximation IS cLS cL high tem- 
perature (small /3) approximation. In any case, this approximation allows me to do 
calculations outside the framework of perturbation theory, but still in leading order in 
a. That is to say, all correction to the effective theory are order of a? or higher (these 
corrections may for instance induce corrections of order a^loga to AT,,). 

Q: In the ratio of IST^jTc oc n^/^, both sides are dimensionless. But A depends on T. If 
you take the ratio, A disappears or not? 

A: If you go through the calculation you will see that the factor n^l^ can be read as 1/A. And 
indeed A depends on T. But the temperature is fixed in the calculation of ^Tc/Tc. 
That is, A is here Ac which can be estimated by taking the critical temperature in 
the absence of interaction. So this temperature dependence does not spoil the linear 
relation between ^Tc/Tc and a. 

Q: OK. Thank you. 
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FIG. 39: Self-energy in the large A'" expansion 

Let me tell you a bit more before ending this lecture. We have established the relation 

AT, 



c 



in which c is a positive constant, given by an integral of the form 



(an^/=^) , (274) 



oo 



coc / dx ^2 ( V (275) 
Jo + cr[x) 

I remind you that a{x) is essentially Sc;(p) — Scj(O). So, provided you know how to calculate 
the self-energy and its momentum dependence, then in principle we can calculate c. The 
problem is to calculate (j{x). 

To do so, we must use non perturbative techniques, of which there are not so many. The 
first technique that I want to discuss is the large N approximation. The idea is to replace the 
field ipo{r), which is a complex field written as 0i(r) +202('") = (0i, 02), by an A^-component 
vector (f){r) = (01, 4>2, - ■ ■ , (Pn)- Then we can write an action 



S= £r{ 0(r) 0(r) + ^ 0(r) , (276) 




which is invariant under 0{N) transformations. The usefulness of this strategy is that 
an analytic calculation is possible when — oo. For the calculation of the self-energy, 
this amounts essentially to resum the chain of bubbles in Fig. 39 (this is actually a. 1/N 
correction; the leading order correction is the mean field correction that leads to a momentum 
independent self-energy). If you remember what I said in the last lecture, this chain of 
bubbles is what produces screening of the long wavelength density fluctuations. Because of 
this screening, the resulting calculation of the self-energy is infrared finite. The calculation 
of the self-energy in this order can be done analytically, and yields c = 2.3. 

Another possible strategy is to do a "brute force" lattice calculation, as you do in QCD 
for instance. This leads to the most accurate determination of c. The results obtained 
almost simultaneously by two groups are c = 1.32 ± 0.02 [18] , and c = 1.29 ± 0.05 [19]. 
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1 X 

(p ~ kc) 

FIG. 40: Schematic behavior of the integrand in Eq. (271) 

After these results were obtained, the calculation of c has become a play ground for 
testing various approximation in field theory or many body physics (see for instance [14] for 
references) . 

Perhaps before closing I should tell you one word about why it is hard to calculate a{x). 
The reason can be understood from looking at the plot of the integrand in Eq. (271). As can 
be seen on Fig. ??, the integrand is peaked around x ~ 1, that is at the momentum p ~ ~ 
a/A^, which delineates the frontier between two very distinct regimes: At momenta smaller 
than Kc, we have the critical regime, well described by the theory of critical phenomena 
which allows in particular a precise determination of the anomalous dimension rj. This is 
well under control. The other regime, that of large momenta, is also very well under control 
because in this regime perturbation theory can be applied (perturbation theory is accurate 
at large momenta in a three dimensional scalar field theory). But the quantity which we 
need is sensitive to what happens at the border line between these two regimes. What we 
need to get with accuracy is the precise point where this transition between the two regimes 
occurs. And this is hard. That is why we need non-trivial non-perturbative techniques to 
do the explicit calculation (and why simple methods, like the large expansion, are off 
quantitatively by almost a factor two). 

What was going to come after that was a discussion of this problem, from the perspective 
of the non-perturbative (or "exact") renormalization group [14]. But that will be for another 
occasion... as the time for these lectures is over. 
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